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Therapeutic  Mechanism  of  BET  Bromodomain  Inhibitor  in  Breast  Cancer 
Abstract 

Genomic  profiling  of  the  MCF7  breast  cancer  cell  line  using  BRD4  ChIP-seq  and  DNase- 
seq  revealed  BRD4  ChIP-seq  and  DNase-seq  to  be  highly  similar.  We  discovered  that 
the  BET  inhibitor  JQ1  tends  to  be  more  effective  in  slowing  the  growth  of  basal  rather 
than  luminal  breast  cancer  cell  lines.  The  initial  gene  expression  response  of  a  basal 
breast  cancer  cell  line,  SUM159,  on  treatment  with  JQ1,  is  predominated  by  the  down- 
regulation  of  gene  expression  and  this  down-regulation  is  highly  associated  with  BRD4 
occupied  genomic  loci.  By  24h  of  JQ1  treatment  secondary  effects  dominate  the  gene 
expression  pattern  and  BRD4  binding  no  longer  predicts  gene  expression  changes. 
Long-term  treatment  of  SUM159  with  JQ1  results  in  this  cell  line  becoming  insensitive  to 
JQ1  and  adopting  a  gene  expression  pattern  that  is  more  luminal-like.  Long-term 
exposure  of  SUM159  cells  to  JQ1  leads  to  drug  resistance  that  is  acquired  through  a 
mechanism  that  allows  BRD4  to  bind  to  chromatin  in  a  bromodomain  independent 
manner. 

Introduction 

Epigenetic  drugs  that  interact  with  chromatin  regulators  to  modify  the  chromatin  states  of 
cancer  cells  promise  to  be  important  in  cancer  treatment  (Yoo  and  Jones,  2006;  Baylin 
and  Jones,  2011).  To  effectively  use  these  drugs  to  manipulate  transcription  for 
therapeutic  purposes,  we  need  to  gain  insight  into  their  genome-wide  function.  We 
proposed  to  develop  both  experimental  and  computational  methods  and  apply  them  to 
identify  the  efficacy,  mechanism,  and  targets  of  BET  domain  inhibitor  JQ1  in  different 
breast  cancer  cells.  This  could  potentially  inform  personalized  breast  cancer  treatment 
based  on  tumor  subtypes.  The  proposed  experimental  and  computational  approaches 
will  be  widely  applicable  to  the  development  of  other  epigenetic  drugs  being  developed 
to  treat  breast  cancers. 

Aim  1 :  Identify  the  cistrome,  transcriptome  and  chromatin  dynamics  in  a  breast  cancer 
cell  line  SUM159  treated  with  BET  domain  inhibitor  JQ1 . 

Aim  2:  Computationally  infer  and  experimentally  validate  the  JQ1  induced  transcriptional 
regulatory  network  and  pathway  in  SUM  159. 

Aim  3:  Apply  the  above  approaches  to  test  JQ1  on  different  breast  cancer  cell  lines  and 
xenografts  to  compare  the  efficacy  and  specificity  of  JQ1  in  different  breast  cancer  cells. 

Keywords 

Breast  cancer,  bromodomain  and  extra  terminal  domain  inhibitors,  JQ1,  triple  negative 
breast  cancer,  luminal  breast  cancer,  basal  breast  cancer,  BRD4,  ChIP-seq 

Accomplishments 

Task  1.  Identify  cistrome,  transcriptome  and  chromatin  dynamics  in  the  breast  cancer 
cell  line  SUM159  treated  with  BET  domain  inhibitor  JQ1  or  BRD4  knockdown. 

To  optimize  the  DNase-seq  protocol,  we  digested  2.5x1 06  breast  cancer  MCF-7  cells 
with  50  units  of  DNase  I  (Roche),  and  then  selected  DNA  fragments  in  50-100,  100-200, 
and  200-300bp  sizes  from  each  of  the  digested  samples  and  barcoded  to  multiplex 
HiSeq  sequencing.  We  carried  out  an  analysis  comparing  the  ability  of  data  sets  from 
different  size  fractions  to  discover  transcription  factor  binding  sites  that  had  been 
identified  by  ChIP-seq.  We  found  that  for  the  estrogen  receptor  and  CTCF  short 
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fragments  were  optimal  (Fig.  1).  Our  analysis  of  the  location  of  DNasel  cleavage  sites 
with  respect  to  nucleosome  positions  revealed  that  chromatin  structure  is  likely  to  be  the 
reason  for  the  effectiveness  of  the  short  fragment  DNase-seq. 


Figure  1.  Fraction  of  estrogen  receptor  (ER)  and  CTCF  ChIP-seq  peaks  in  MCF7 
identified  as  DNase-seq  peaks  for  different  fragment  size  selections.  Short  fragments 
are  much  more  efficient  than  longer  ones  in  recovering  CTCF  binding  sites  but  this  effect 
is  less  pronounced  for  the  estrogen  receptor. 

We  determined  that  DNase-seq  data  looks  highly  similar  to  BRD4  ChIP-seq  and  did  not 
provide  much  new  information  about  the  influence  of  JQ1  on  chromatin  and 
transcriptional  changes  (Fig  2).  This  is  a  consistent  with  the  known  role  of  BRD4  in 
opening  chromatin  (Voigt  and  Reinberg,  2011),  although  the  degree  of  similarity 
between  BRD4  ChIP-seq  and  DNase-seq  is  remarkable.  In  subsequent  analyses  we  will 
use  H3K27ac,  a  mark  that  is  indicative  of  enhancer  activity,  in  addition  to  BRD4  ChIP- 
seq. 
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Figure  2.  BRD4  ChIP-seq  and  DNase-seq  in  untreated  MCF-7  cells  and  after  1h  and  6h 
of  JQ1  treatments.  BRD4  ChIP-seq  signal  is  very  similar  to  DNasel  hypersensitivity  and 
does  not  provide  much  complementary  information. 


Task  2:  Computationally  infer  and  experimentally  validate  the  JQ1  induced 
transcriptional  regulatory  network  and  pathway  in  SUM159. 

We  conducted  growth  rate  analyses  on  a  panel  of  breast  cancer  cell  lines  treated  by 
different  concentrations  of  JQ1 .  This  analysis  revealed  that  basal  breast  cancer  cell  lines 
tend  to  be  more  responsive  to  JQ1  treatment  than  luminal  cell  lines  (Fig.  3).  For  this 
reason  we  decided  to  focus  on  the  basal  SUM159  cell  line  instead  of  the  luminal  MCF-7 
cell  line. 
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Very  sensitive  Resistant 


Figure  3:  Response  to  BET  bromodomain  inhibition  in  breast  cancer.  Heatmap  of 
mean  IC50s  of  BET  inhibitors  and  inactive  analogues  (JQ1R  and  alprazolam)  as 
negative  controls  in  a  panel  of  TNBC  (red),  HER2+  (light  blue),  and  luminal  (blue)  breast 
cancer,  and  immortalized  basal  mammary  epithelial  (orange)  cell  lines.  IC50  values  are 
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represented  in  red-white-black  scale:  very  sensitive  (<0.01  pM)  in  red,  sensitive  (=1pM) 
in  white,  to  resistant  (>20pM)  in  black. 

We  analyzed  the  genome-wide  gene  expression  profile  of  SUM159  after  JQ1  treatment 
at  3h,  12h  and  24h  in  relation  to  a  control  group  treated  with  DMSO.  This  analysis 
revealed  that  genes  were  predominantly  down-regulated  after  short  term  (3h)  treatment 
by  JQ1,  as  expected  from  current  models  of  JQ1  action  (Filippakopoulos  et  al.,  2010). 
We  used  a  model  of  “regulatory  potential”  to  use  BRD4  to  predict  the  genes  likely  to  be 
regulated  by  BRD4  and  effected  by  JQ1  treatment.  In  this  model  a  score  is  calculated  for 
each  gene  as  a  sum  of  weights  in  which  BRD4  ChIP-seq  peaks  close  to  the  transcription 
start  site  of  the  gene  are  given  more  weight  than  BRD4  peaks  that  are  far  from  the  TSS 
of  that  gene.  We  evaluated  the  performance  of  this  model  in  predicting  gene  expression 
changes  associated  with  3h,12h  and  24h  of  JQ1  treatment  and  found  that  BRD4  ChIP- 
seq  could  inform  which  genes  were  down-regulated  at  3h  but  could  not  predict  the  up- 
regulated  genes  (Fig  4).  The  performance  of  this  model  deteriorated  with  treatment  time 
as  secondary  effects  dominated  the  primary  effect  of  BET  inhibition.  A  similar  analysis 
based  on  H3K27ac  ChIP-seq  data  resulted  in  a  poorer  prediction  for  down -regulated 
genes  but  also  some  predictive  power  for  up-regulated  genes. 


Figure  4.  Regulatory  potential  model  performance  using  BRD4  ChIP-seq  to  predict 
differentially  expressed  genes  for  3h  (A),  12h  (B),  and  24h  (C)  of  JQ1  treatment.  BRD4 
ChIP-seq  can  predict  early  (3h)  down-regulated  genes  but  not  up  genes  or  late  (24h) 
genes. 

We  analyzed  gene  expression  levels  of  3hr  JQ1  up-,  down-  and  non-regulated  genes  in 
a  large  collection  of  luminal  A,  luminal  B  and  basal  breast  cancer  gene  expression 
profiles  (METABRIC)  (Fig  5A).  JQ1  down-regulated  genes  are  more  highly  expressed  in 
basal  cells  than  in  either  of  the  luminal  cell  types.  This  trend  suggests  that  JQ1  down- 
regulates  the  genes  that  are  associated  with  the  basal  nature  of  SUM159  pushing  the 
gene  expression  pattern  towards  the  luminal  type.  After  long-term  exposure  to  JQ1 
SUM159  cells  acquire  resistance  to  this  BET  inhibitor.  Differential  gene  expression 
analysis  comparing  SUM159R  with  SUM159  shows  that  genes  that  are  down-regulated 
in  SUM159R  are  more  highly  expressed  in  the  basal  cell  type  (Fig  5B).  This  suggests 
that  the  early  effects  of  BET  inhibition  that  make  SUM159  cells  more  luminal-like  do  not 
disappear  on  acquisition  of  resistance. 
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Figure  5.  BET  bromo- 
domain  inhibition 
gene  expression 
response.  Down-,  up- 
and  non-regulated 
genes  were  identified 
from  RNA-seq  data 
comparing  SUM159 
cells  after  3hr  of  JQ1 
treatment  with  3hr  of 
DMSO.  Boxplots 

represent  the  median 
and  quartiles  of 
standardized  gene 
expression  in 

METABRIC  basal, 
luminal  A  and  luminal  B 


data  (A).  Gene  down- 
regulated  by  JQ1  tend  to  be 
more  highly  expressed  in  basal  cells.  In  (B)  differential  genes  are  identified  from  a 
comparison  of  SUM159  and  BET  resistant  SUM159  (SUM159R)  cells  under  3hr  DMSO. 
The  genes  that  are  down-regulated  are  more  highly  expressed  in  basal  cell. 


To  further  investigate  the  nature  of  the  BET  resistant  cells  grown  under  low 
concentrations  of  JQ1  we  exposed  them  to  a  pulse  dose  of  higher  concentration  JQ1 
and  analyzed  the  gene  expression  response  relative  to  cells  treated  with  DMSO  at  3hr, 
12hr  and  24hr.  The  direction  of  response  of  differentially  expressed  genes  was  highly 
consistent  between  SUM159  and  SUM159R  cells  at  all  time  points  (Fig  6).  This  shows 
that  gene  expression  in  SUM159R  cells  remains  under  some  control  of  bromodomain 
associated  Brd4. 
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Figure  6.  Gene  expression  analysis  by  mRNA-seq  in  SUM159  and  JQ1  resistant 
SUM159  (SUM159R)  cells  3hr,  12hr  and  24hr  after  JQ1  treatment.  Differentially 
expressed  genes  associated  with  JQ1  treatment  in  both  SUM159  and  SUM159R  are 
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shown  in  red.  The  direction  of  response  of  genes  is  consistent  between  SUM159  and 
SUM159R  cells. 


We  used  immunofluorescence  to  examine  the  basal/luminal  nature  of  SUM159  and 
MDA-MB-231  cells  treated  for  3  days  by  JQ1  or  DMSO  as  a  control.  These  data  are 
consistent  with  the  gene  expression  analysis,  showing  the  acquisition  of  luminal 
cytokeratin  markers  and  the  loss  of  basal  ones  on  JQ1  treatment  (Fig  7). 


SUM159  MDA-MB-231 

DMSO  JQ1  DMSO  JQ1 


Figure  7.  Immunofluorescence  of 
basal  (high  molecular  weight  -HMW) 
and  luminal  (low  molecular  weight  - 
LMW  and  cytokeratin  18) 
cytokeratins  in  SUM159  and  MDA- 
MB-231  TNBC  lines  before  and  3 
days  after  JQ1  (500nM)  treatment. 
Scale  bars  show  20pm. 


Task  3:  Apply  the  above  approaches  to  test  JQ1  on  different  breast  cancer  cell  lines  and 
xenografts  to  compare  the  efficacy  and  specificity  of  JQ1  in  different  breast  cancer  cells. 


We  have  collaborated  with  Dr  Polyak’s  lab  to  assess  the  efficacy  of  JQ1  on  different 
breast  cancer  cells  lines  and  discovered  that  JQ1  tends  to  slow  the  growth  of  basal  cell 
types  more  than  luminal  ones.  In  this  collaboration  we  have  also  derived  a  JQ1 
insensitive  cell  line  (SUM159R)  derived  from  SUM159  cells  subjected  to  long-term 
treatment  with  JQ1.  We  show  using  10  independent  single  cell-derived  clones  that  this 
insensitive  cell  line  is  not  likely  due  to  the  outgrowth  of  a  subpopulation  of  pre-existing 
resistant  cells  (Fig.  8A).  Xenografts  derived  from  SUM159R  cells  were  not  responsive 
to  JQ1  treatment  in  parallel  studies  with  sensitive  TNBC  cells  (Fig.  8B). 
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Figure  8  A,  Cellular  viability  of  SUM159,  and  pool  and  single  cell  clones  of  SUM159R 
cells  treated  with  different  doses  of  JQ1.  B,  Tumor  weight  of  xenografts  derived  from 
SUM159  and  SUM159R  cells.  Mice  were  administered  JQ1  for  14  (SUM159)  and  30 
(SUM159R)  days  beginning  at  day  14  and  26,  respectively,  after  injection. 


Genome-wise  CRISPR-Cas9  screens  of  genes  in  parental  and  resistant  SUM159  cell 
lines  are  underway.  For  these  data  we  will  assess  which  genes  support  or  hinder  cell 
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proliferation  in  the  presence  and  absence  of  JQ1 .  Based  on  the  results  of  the  screen  we 
will  assess  the  effect  of  gene  knock-outs  on  genome-wide  gene  expression  to  examine 
the  extent  to  which  these  knock-outs  can  modulated  the  basal  to  luminal  transition  and 
thereby  yield  insight  into  the  higher  sensitivity  of  basal  breast  cancer  cell  types  to  the 
BET  inhibitor  JQ1 . 

Key  Research  Accomplishments 

•  Optimized  DNase-seq  conditions  and  developed  a  model  that  explains  fragment 
size  effects  in  DNase-seq  data. 

•  Analyzed  and  integrated  BRD4  ChIP-chip  and  expression  data  in  SUM159  cell 
line,  before  and  after  JQ1  treatment. 

•  Conducted  H3K27ac  histone  mark  ChIP-seq  experiments  in  SUM159  cell  lines. 

•  Developed  a  methodology  for  the  systematic  analysis  of  relationships  between 
H3K27ac  profiles  and  gene  expression. 

•  Established  a  JQ1  insensitive  cell  line  from  a  parental  cell  line  that  is  JQ1 
sensitive. 

•  Identified  factors  that  underlie  JQ1  resistance. 

•  Manuscript:  “Response  and  resistance  to  BET  bromodomain  inhibitors  in  triple 
negative  breast  cancers.”  Shu  et  al,  accepted  by  Nature,  2016 

Impact 

Here  we  provide  support  for  BET  inhibition  in  the  treatment  of  triple  negative  breast 
cancer.  We  discovered  resistance  to  BET  inhibition  in  which  BRD4  is  recruited  to 
chromatin  in  a  bromodomain  independent  manner.  These  findings  suggest  that  BET 
inhibition  may  need  to  be  combined  with  complementary  drugs  to  address  mechanisms 
of  resistance. 

Changes/Problems 

We  originally  proposed  a  study  of  BET  inhibition  in  the  luminal  breast  cancer  cell  line 
MCF7.  We  discovered  that  luminal  breast  cancers  tend  to  be  less  sensitive  to  BET 
inhibition  than  triple  negative  ones.  We  turned  our  focus  to  triple  negative  breast 
cancers  as  these  do  not  currently  have  good  drug  treatment  options  and  BET  inhibition 
appears  to  be  a  promising  route. 

Products 

Manuscript:  “Response  and  resistance  to  BET  bromodomain  inhibitors  in  triple  negative 
breast  cancers.”  Shu  et  al,  accepted  by  Nature,  2016 
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Response  and  resistance  to  BET  bromodomain 
inhibitors  in  triple -negative  breast  cancer 

Shaokun  Shu1,2*,  Charles  Y.  Lin1,2*,  Housheng  Hansen  He1,2,3,4,5*,  Robert  M.  Witwicki1,2*,  Doris  R  Tabassum1,  Justin  M.  Roberts1, 
Michalina  Janiszewska1,2,  Sung  Jin  Huh1,2,  Yi  Liang4,  Jeremy  Ryan1’2,  Ernest  Doherty1,6,  Hisham  Mohammed7,  Hao  Guo3, 

Daniel  G.  Stover1,2,  Muhammad  B.  Ekram1,2,  Guillermo  Peluffo1,2,  Jonathan  Brown1,2,  Clive  D’Santos7,  Ian  E.  Krop1,2, 

Deborah  Dillon1,8,  Michael  McKeown1,2,  Christopher  Ott1,2,  Jun  Qi1,2,  Min  Ni1,2,  Prakash  K.  Rao9,  Melissa  Duarte9, 

Shwu-Yuan  Wu10,  Cheng-Ming  Chiang10,  Lars  Anders11,  Richard  A.  Young11,  Eric  P.  Winer1,2,  Antony  Letai1,2,  William  T.  Barry2,3, 
Jason  S.  Carroll7,  Henry  W.  Long1,9,  Myles  Brown1,2,9,  X.  Shirley  Liu3,9,12,  Clifford  A.  Meyer3,  James  E.  Bradner1,2,12  & 

Kornelia  Polyak1,2,9,12 


Triple-negative  breast  cancer  (TNBC)  is  a  heterogeneous  and 
clinically  aggressive  disease  for  which  there  is  no  targeted  therapy1-3. 
BET  bromodomain  inhibitors,  which  have  shown  efficacy  in  several 
models  of  cancer4-6,  have  not  been  evaluated  in  TNBC.  These 
inhibitors  displace  BET  bromodomain  proteins  such  as  BRD4 
from  chromatin  by  competing  with  their  acetyl-lysine  recognition 
modules,  leading  to  inhibition  of  oncogenic  transcriptional 
programs7-9.  Here  we  report  the  preferential  sensitivity  of  TNBCs 
to  BET  bromodomain  inhibition  in  vitro  and  in  vivo ,  establishing 
a  rationale  for  clinical  investigation  and  further  motivation  to 
understand  mechanisms  of  resistance.  In  paired  cell  lines  selected 
for  acquired  resistance  to  BET  inhibition  from  previously  sensitive 
TNBCs,  we  failed  to  identify  gatekeeper  mutations,  new  driver 
events  or  drug  pump  activation.  BET-resistant  TNBC  cells  remain 
dependent  on  wild-type  BRD4,  which  supports  transcription 
and  cell  proliferation  in  a  bromodomain-independent  manner. 
Proteomic  studies  of  resistant  TNBC  identify  strong  association 
with  MED1  and  hyper-phosphorylation  of  BRD4  attributable  to 
decreased  activity  of  PP2A,  identified  here  as  a  principal  BRD4 
serine  phosphatase.  Together,  these  studies  provide  a  rationale 
for  BET  inhibition  in  TNBC  and  present  mechanism-based 
combination  strategies  to  anticipate  clinical  drug  resistance. 

To  explore  non-oncogene  addiction  to  BRD4  in  breast  cancer,  we 
studied  a  series  of  BET  bromodomain  inhibitors  (BBI)  across  breast 
cell  lines  reflecting  transcriptionally  defined  breast  cancer  subtypes: 
luminal,  HER2+  and  TNBC2,10,  as  well  as  MCF10A  and  MCF12A  basal/ 
mesenchymal  immortalized  mammary  epithelial  cells  (Supplementary 
Table  1).  Potent  inhibitory  effects  were  observed  preferentially  in 
TNBC  lines,  compared  to  more  resistant  luminal  lines  (Fig.  la). 
Analysis  of  potency  of  drug  response  and  subtype  or  known  driver 
mutations  identified  the  basal  subtype  as  the  only  significant  associa¬ 
tion  (P= 0.0475)  (Supplementary  Table  1  and  data  not  shown).  BRD4 
dependency  was  confirmed  by  RNA  interference  and  phenocopied  BBI 
(Extended  Data  Fig.  la-c).  JQ1  or  BRD4  knockdown  induced  growth 
inhibition  and  resulted  in  G1  arrest  and  apoptosis  (Extended  Data 
Fig.  ld-g).  Expression  of  factors  described  to  mediate  JQ1  effect  (MYC) 
or  required  for  TNBC  growth  (JAK2/STAT3)  showed  no  clear  associ¬ 
ation  with  JQ1  sensitivity  (Extended  Data  Fig.  lh  and  Extended  Data 
Fig.  2a,  b).  JQ1  treatment  of  TNBC  cells  induced  significant  morphologic 


changes  consistent  with  induction  of  senescence,  confirmed  by 
(3-galactosidase  staining  and  luminal  differentiation  evidenced  by 
changes  in  the  expression  of  basal  and  luminal  markers  (Extended 
Data  Fig.  2c,  d  and  Fig.  lb). 

Extending  the  translational  significance  of  these  findings,  we  evalu¬ 
ated  the  ability  of  JQ1  to  inhibit  tumour  growth  in  murine  TNBC  xen¬ 
ografts.  Two-week  treatment  efficiently  inhibited  established  tumour 
growth  from  SUM  159  and  MDA-MB-231  lines,  and  patient-derived 
primary  human  TNBC  xenografts  (Fig.  lc  and  Extended  Data  Fig.  2e,  f). 
Downregulation  of  BRD4  using  two  independent  TET-inducible  short 
hairpin  RNAs  produced  even  more  pronounced  effects,  leading  to  com¬ 
plete  tumour  regression  and  failure  to  regrow  even  after  discontinuing 
doxycycline  treatment  (Fig.  lc  and  Extended  Data  Fig.  2g).  Evidence 
of  BBI-inducedbasal-to-luminal  differentiation  was  confirmed  in  vivo 
(Extended  Data  Fig.  2f,  h). 

Using  integrated  epigenomic  analysis  (Supplementary  Table  2), 
we  identified  the  direct  transcriptional  targets  of  BBI  in  TNBC.  BBI 
binding  was  identified  at  active  promoter  and  enhancer  regions 
using  ChemSeq11  for  biotinylated  JQ1  (Bio-JQl)  enrichment  and 
chromatin  immunoprecipitation  followed  by  sequencing  (ChIP-seq) 
for  acetyl-histone  (H3K27ac)  and  BRD4  enrichment,  with  the  three 
marks  showing  near  perfect  co -localization  (Fig.  Id  and  Extended  Data 
Fig.  3a).  BBI  efficiently  displaced  chromatin-bound  BRD4  in  treated 
SUM  159  (Fig.  le  and  Extended  Data  Fig.  3b)  and  in  SUM  149  cells 
(Extended  Data  Fig.  3c).  To  identify  biologically  relevant,  direct  tar¬ 
gets  of  BBI  in  SUM  159  and  SUM  149  cells,  we  quantified  binding  of 
Bio-JQl  and  BRD4  genome-wide  and  found  strong  enrichment  at  219 
and  159  super-enhancers,  respectively  (Fig.  If,  Extended  Data  Fig.  3d 
and  Supplementary  Table  3)8>9>12>13.  Transcription  factors  with  known 
roles  in  breast  cancer,  such  as  POU5F1B/MYC14  and  HIFla15,  were  evi¬ 
dent  among  top  super-enhancer-associated  genes  in  both  lines.  Kinetic 
effects  of  JQ1  treatment  on  gene  expression  demonstrated  preferential 
super-enhancer- associated  gene  downregulation  (Fig.  lg  and  Extended 
Data  Fig.  3e,  f).  Expression  changes  were  observed  within  3  h  after  JQ1 
treatment  and,  as  expected,  more  genes  were  significantly  down-  than 
upregulated  (Extended  Data  Fig.  3g-j  and  Supplementary  Table  4). 
Unsupervised  Metacore16  analysis  of  JQ1- affected  target  pathways 
revealed  downregulation  of  regulatory  and  effector  genes  in  anti- 
apoptotic  and  JAK/STAT  signalling  pathways  (Extended  Data  Fig.  3k). 
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Figure  1  |  Response  to  BBIs  in  breast  cancer,  a,  Heat  map  of  mean 
IC50s  of  BBIs  and  inactive  analogues  in  breast  cell  lines.  Error  bars 
represent  s.e.m.  b,  Immunofluorescence  of  basal  (basCK)  and  luminal 
(lumCK  and  CK18)  cytokeratins  in  TNBC  lines.  Scale  bars,  20  pm. 
c,  Box  plots  depicting  xenograft  weights;  n  indicates  the  number  of  mice 
per  experiment.  *P<  0.0001  (unpaired  f-test).  d,  Heat  map  showing 
biotinylated  JQ1  (Bio-JQI),  BRD4  and  H3K27ac  binding  at  transcription 
start  site  (TSS)  and  Bio-JQI -bound  enhancer  regions.  Each  row  represents 
a  single  genomic  region  (±10  kb)  from  TSS  or  enhancer  centre.  Genomic 
occupancy  is  shaded  by  binding  intensity  in  units  of  reads  per  million  per 
base  pair  (rpm  per  bp),  e,  Gene  tracks  depicting  Bio-JQI  and  BRD4  with 
or  without  JQ1  in  SUM  159  cells  at  the  HIF1A  locus,  x-axis,  chromosome 
position  with  gene  structures  below;  y- axis,  genomic  occupancy  in  units  of 
rpm  per  bp;  red  bar,  HIF1A  super-enhancer,  f,  Plot  of  enhancers  defined  in 
untreated  SUM159  cells  ranked  by  increasing  Bio-JQI  signal  (units  rpm). 
Grey  line  marks  cut-off  discriminating  typical  from  super-enhancers, 
g,  Box  plots  showing  the  log2  fold  change  in  expression  relative  to  control 
of  either  all  active  or  super- enhancer  (SE) -associated  genes  upon  JQ1 
treatment.  **P<  10-5  (Welsh’s  f-test);  *P<  10-3  (Welsh’s  f-test). 


These  data  support  selective  disruption  of  super-enhancer-associated 
genes  by  JQ1,  leading  to  deregulation  of  coordinated  transcriptional 
pathways  involved  in  cell  proliferation,  invasion  and  survival. 

Dissecting  resistance  to  targeted  therapy  is  critical  to  elucidate  mech¬ 
anisms  of  drug  and  target  action,  and  to  suggest  approaches  to  treat 
or  anticipate  drug  resistance  in  patients.  Therefore,  we  established 
BBI-resistant  TNBC  cell  lines  by  long-term  culture  of  both  SUM  159 
and  SUM149  cells  in  escalating  JQ1  doses.  Low  (0.5  pM)  and  high 
(2.0  pM)  doses  of  JQ1  severely  impaired  proliferation  of  parental 
SUM  159  and  SUM  149  lines,  reducing  viable  cells  after  6  days  (Fig.  2a 
and  Extended  Data  Fig.  31).  In  contrast,  JQl-resistant  cells  (SUM159R 
and  SUM149R)  proliferated  linearly,  even  in  high  JQ1  doses  (20  pM) 
(Fig.  2a  and  Extended  Data  Fig.  31).  BBI-resistance  is  not  attributable 
to  drug  export,  as  MDR1  and  other  transporters  are  not  transcription¬ 
ally  upregulated  (Extended  Data  Fig.  4a),  co-incubation  with  MDR1 
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Figure  2  |  Acquired  BBI-resistance  in  TNBC.  All  error  bars  represent 
s.d.,  n  =  3.  a,  Viable  cell  numbers  after  JQ1  treatment.  *P  <  0.0001  (two- 
way  ANOVA).  b,  Cellular  viability  after  treatment  with  BBIs.  *P  <  0.0001 
(non-linear  regression,  extra  sum-of- squares  test),  c,  Genomic  regions 
containing  a  super-enhancer  in  SUM159  or  SUM159R  cells  ranked  by 
log2  change  in  Bio-JQI  genomic  binding  signal,  x-axis,  log2  fold  change  in 
Bio-JQI  signal  coloured  by  intensity  of  change,  d,  e,  Box  plot  showing  the 
log2  fold  change  in  BRD4  genomic  occupancy  (d)  and  gene  expression  (e) 
at  regions  with  gained,  conserved,  or  lost  Bio-JQI  binding  in  SUM159R 
versus  SUM159  cells,  f,  g,  Gene  tracks  depicting  Bio-JQI,  BRD4,  and 
H3K27ac  at  the  BCL-xL  (also  known  as  BCL2L1;  f)  and  SOD2  (g)  locus. 
The  x-axis  shows  position  along  the  chromosome  with  gene  structures 
drawn  below.  They-axis  shows  genomic  occupancy  in  units  of  rpm  per  bp. 
h,  Box  plot  showing  the  log2  fold  change  in  BRD4  genomic  occupancy  at 
regions  bound  by  Bio-JQI. 


inhibitors  (verapamil)  had  no  effect  (Extended  Data  Fig.  4b),  and  struc¬ 
turally  divergent  BBIs  are  equally  inactive  as  JQ1  (Fig.  2b).  Further 
support  is  provided  by  the  equivalent  chromatin  engagement  of  BRD4 
in  sensitive  and  resistant  cells,  demonstrated  by  binding  with  Bio-JQI 
(Extended  Data  Fig.  4c).  Notably,  BBI-resistant  TNBC  cells  retain 
sensitivity  to  compounds  from  orthogonal  active  drug  classes,  such  as 
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CXCR2  and  JAK2  inhibitors17,  establishing  specific  resistance  to  BBIs 
(Extended  Data  Fig.  4d).  Adaptive  drug  resistance  was  not  attributable 
to  outgrowth  of  a  minor  subpopulation  of  pre-existing  resistant  cells, 
as  10  independent  single-cell-derived  clones  showed  similar  resistance 
profiles  to  pooled  SUM159R  cells  (Extended  Data  Fig.  4e).  Similar 
results  were  obtained  in  vivo ,  as  SUM159R-derived  xenografts  were 
JQ1 -unresponsive  (Extended  Data  Fig.  4f).  In  all  resistant  TNBC  pop¬ 
ulations  studied,  exome  sequencing  failed  to  identify  alterations  in  BET 
bromodomain-encoding  genes  (for  example,  gatekeepers)  or  known 
driver  genes  (parallel  pathway  activation)  (Supplementary  Table  5). 

Absent  new  genetic  alterations,  we  explored  the  plausibility  of  an 
epigenomic  mechanism  of  resistance.  Differential  enhancer  analysis 
revealed  a  significant  gain  of  super-enhancers  in  resistant  SUM159R 
cells  (Fig.  2c  and  Supplementary  Table  6).  The  gain  of  Bio-JQl  super 
enhancers  was  associated  with  enrichment  for  BRD4  binding  to  these 
genomic  loci  (Fig.  2d)  and  increased  transcription  of  associated  genes 
(Fig.  2e).  An  upstream/intragenic  region  of  H3K27ac  at  the  BCL-xL 
locus  featured  prominently  among  top  gained  super  enhancers  in 
SUM159R  (Fig.  2f),  consistent  with  increased  BCF-xF  messenger 
RNA  and  protein  expression  in  resistant  cells  (Supplementary  Table  7, 
Extended  Data  Fig.  4g).  Functionally,  cells  with  acquired  resistance  to 
BBI  featured  a  concordant  switch  in  JQ1  anti-apoptotic  response  based 
on  dynamic  BH3  profiling18,19  (Extended  Data  Fig.  4h). 

Observing  emergent  enhancers  in  resistant  cells,  we  assessed  whether 
BBI-resistant  TNBC  cells  retained  non-oncogene  addiction  to  BRD4. 
Notably,  we  observed  loss  of  SUM159R  cell  viability  upon  BRD4 
knockdown  (Extended  Data  Fig.  5a,  b).  Together  these  studies  estab¬ 
lish  persistence  of  BRD4  addiction  despite  resistance  to  bromodomain 
inhibition,  establishing  the  plausibility  of  bromodomain-independent 
recruitment  of  BRD4  to  enhancers  in  BBI-resistant  TNBCs.  To  test  this 
hypothesis,  we  performed  BRD4  ChIP-seq  on  sensitive  and  resistant 
cells  with  and  without  JQ1.  JQ1  neither  displaced  BRD4  from  chro¬ 
matin  in  SUM159R  (Fig.  2g),  nor  meaningfully  influenced  epigenome 
structure  by  H3K27ac  ChIP-seq  (Extended  Data  Fig.  5c-g).  Notably, 
several  luminal  markers  (FOXA1,  CD24,  and  luminal  cytokeratins) 
were  elevated  in  SUM159R  cells  in  cell  culture  and  in  vivo  (Extended 
Data  Fig.  5h,  i),  supporting  a  model  whereby  resistance  arises  via  essen¬ 
tial  BRD4  recruitment  to  chromatin  in  a  bromodomain-independent 
manner.  Similar  observations  were  made  in  SUM149R  cells  and  in 
TNBC  cells  inherently  resistant  to  JQ1  (Extended  Data  Fig.  3h-j; 
Extended  Data  Fig.  6a-d),  suggesting  a  general  mechanism  of  epig¬ 
enomic  resistance  to  BBI. 

To  disclose  potential  differences  in  BRD4-associated  complexes 
between  sensitive  and  resistant  SUM  159  cells,  we  performed  quanti¬ 
tative  proteomics  using  RIME  (rapid  immunoprecipitation  mass  spec¬ 
trometry  of  endogenous  proteins)20  with  and  without  JQ1.  Analysis 
of  BRD4-associated  proteins  identified  relative  enrichment  of  MED  1 
and  BRD3  in  JQ1 -treated  resistant  cells  (Fig.  3a,  Extended  Data  Fig.  7 
and  Supplementary  Table  8).  BRD4  immunoprecipitation  followed 
by  immunoblot  for  MED1  and  BRD3  revealed  that  JQ1  efficiently 
displaced  BRD4  from  MED1  in  sensitive  cells,  but  not  in  resistant  cells 
(Fig.  3b),  a  result  confirmed  in  SUM149  and  BBI-resistant  SUM149R, 
as  well  as  inherently  resistant  TNBC  and  luminal  lines  (Extended 
Data  Fig.  8a).  Although  elevated  BRD3  abundance  was  observed  in 
SUM159R  cells,  increased  association  of  BRD4  and  BRD3  was  not 
confirmed  by  immunoblot  (Fig.  3b).  To  assess  functionally  whether 
increased  recruitment  of  BRD4  to  chromatin  by  MED1  under¬ 
lies  resistance  to  JQ1,  we  expressed  an  exogenous  bromodomain- 
inactivated  mutant  (BDmut)  with  concomitant  knockdown  of 
endogenous  BRD4  (Extended  Data  Fig.  8b,  c).  Downregulation  of 
endogenous  BRD4  decreased  cell  growth  both  in  parental  and  resist¬ 
ant  cells,  which  was  rescued  by  enforced  expression  of  wild-type 
BRD4  (Fig.  3c).  BDmut  BRD4  expression  failed  to  rescue  parental 
SUM  159  cells,  but  supported  growth  of  JQ1 -resistant  SUM159R  con¬ 
sistent  with  an  evident  bromodomain-independent  mechanism  of 
BRD4  recruitment  (Fig.  3c).  Next,  we  assessed  the  sensitivity  of  cells 
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Figure  3  |  Mechanism  of  BBI-resistance  in  TNBCs.  All  error  bars 
represent  s.d.,  n  =  3.  a,  Plot  depicting  changes  in  BRD4-associated  proteins 
in  SUM159  and  SUM159R  cells  following  JQ1  treatment  based  on  SILAC 
RIME.  The  axes  represent  logio  of  fold  change  (FC).  b,  Immunoblot 
analysis  of  BRD4  immunoprecipitates  and  total  cell  lysates  in  SUM159  and 
SUM159R  cells.  For  gel  source  data,  see  Supplementary  Fig.  1.  c,  Cellular 
viability  of  SUM  159  and  SUM159R  cells  expressing  exogenous  wild- 
type,  BDmut,  7  A  and  7D  mutant  BRD4  with  concomitant  knockdown  of 
endogenous  BRD4.  *P  =  0.04  (paired  f-test).  d,  e,  Sensitivity  of  SUM159 
(d)  and  SUM159R  (e)  cells  expressing  exogenous  wild-type  or  BDmut 
BRD4  to  JQ1  with  concomitant  knockdown  of  endogenous  BRD4. 

*P  <  0.0001  (non-linear  regression,  extra  sum-of-squares  test). 


expressing  BDmut  BRD4  to  JQ1,  observing  increased  sensitivity  to 
JQ1  in  parental  SUM  159  cells  exogenously  expressing  BDmut  (Fig. 
3d).  In  contrast,  expression  of  BDmut  BRD4  in  SUM159R  cells  res¬ 
cued  the  anti -proliferative  effect  of  JQ1  (Fig.  3e),  although  this  could 
partially  be  due  to  the  slower  growth  of  BDmut  expressing  cells. 
Together,  these  studies  suggest  that  BBI-resistance  is  associated  with 
increased  binding  of  BRD4  to  MED1,  in  a  bromodomain-independ¬ 
ent  manner  unaffected  by  JQ1. 

A  recent  study  reported  that  the  stability  and  nuclear  localization  of 
BRD4  is  increased  with  phosphorylation  by  casein  kinase  II  (CK2)21.  To 
explore  the  contribution  of  BRD4  phosphorylation  to  BBI-resistance, 
we  performed  immunoblot  analysis  in  parental  and  resistant  cells 
and  found  a  marked  increase  of  phospho-BRD4  (pBRD4)  in  resistant 
cells  (Fig.  4a  and  Extended  Data  Fig.  8d).  Small-molecule  inhibition 
of  CK2  decreased  BRD4  phosphorylation  in  SUM  159  and  SUM159R 
cells  (Extended  Data  Fig.  8e).  These  results  imply  BRD4  hyperphos¬ 
phorylation  in  resistant  cells  either  owing  to  increased  phosphoryla¬ 
tion  by  CK2  or,  alternatively,  to  decreased  dephosphorylation  by  an  as 
yet  unidentified  BRD4  phosphatase.  We  therefore  first  analysed  CK2 
activity  in  parental  and  resistant  cells  by  performing  pan-CK2  substrate 
immunoblots  and  detected  no  significant  differences  in  CK2  activity 
(Extended  Data  Fig.  8f). 
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Figure  4  |  Regulation  and  relevance  of  BRD4  phosphorylation. 

All  error  bars  represent  s.d.,  n  —  3.  a,  Immunoblot  for  the  indicated 
proteins  following  JQ1  treatment,  b,  Immunoblot  for  the  indicated 
proteins  after  knockdown  of  PP2A  A  or  C  or  both  subunits,  c,  Viable 
cell  numbers  of  JQ1 -treated  control  and  shPP2A-C  expressing  SUM149 
cells.  Left,  *P<  0.0001  (two-way  ANOVA);  right,  *P  =  0.0037and 
**P  =  0.0015  (two-way  ANOVA).  d,  Immunoblot  of  pBRD4  and  BRD4 
in  SUM159R  cells  following  phenothiazine  (PTZ)  treatment,  e,  Viable 
cell  numbers  of  SUM159R  cells  treated  with  JQ1,  phenothiazine  or  both 
compounds.  *P<  0.0001  (two-way  ANOVA).  f,  g,  Immunoblot  of  BRD4 
immunoprecipitates  and  total  cell  lysates  of  SUM159R  cells  after  3  h 
treatment  with  JQ1  and  CK2i  (f)  and  JQ1  and  PTZ  (g).  h,  Immunoblot 
of  Flag-BRD4  (wild  type  or  mutant)  immunoprecipitates  and  total  cell 
lysates  after  3  h  treatment  with  JQ1.  i,  JQ1  sensitivity  of  SUM159  and 
SUM159R  cells  expressing  exogenous  wild-type  or  mutant  (7 A,  7D)  BRD4. 
*P<  0.0001  (non-linear  regression,  extra  sum-of-squares  test).  For  gel 
source  data,  see  Supplementary  Fig.  1. 

Inactivation  of  the  PP2A  phosphatase  tumour  suppressor  gene 
occurs  commonly  in  breast  cancer  and  is  associated  with  therapy 
resistance22;  PP2A  also  often  opposes  CK2  function23,24.  Thus,  we 
investigated  whether  PP2A  may  dephosphorylate  BRD4  and  whether 
decreased  PP2A  activity  could  lead  to  BBI  resistance.  Downregulation 
of  PP2A  catalytic  subunit  ( PPP2CA )  in  SUM  149  and  SUM  159  cells  led 
to  increased  BRD4  phosphorylation,  establishing  PP2A  as  a  previously 
unrecognized  BRD4  phosphatase  (Fig.  4b),  further  supported  by  phar¬ 
macologic  inhibitors  of  PP2 A  that  showed  similar  effects  (Extended 
Data  Fig.  8g).  To  strengthen  the  link  between  PP2A  activity  and  BBI 
resistance,  we  tested  the  JQ1  sensitivity  of  SUM  149  cells  following  the 
knockdown  of  PP2 A  C  subunit  and  determined  that  downregulation 
of  PP2A  decreased  JQ1  sensitivity  (Fig.  4c).  We  have  collaboratively 
reported  phenothiazine  compounds  as  activators  of  PP2 A  enzymatic 
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activity25.  Thus,  we  analysed  pBRD4  levels  in  SUM159R,  SUM149R  and 
other  cell  lines  after  short-term  treatment  with  phenothiazine  (PTZ) 
and  detected  rapid  dephosphorylation  of  BRD4  (Fig.  4d  and  Extended 
Data  Fig.  8h).  Combined  treatment  with  PTZ  and  JQ1  overcame  BBI 
resistance  in  SUM159R  cells  (Fig.  4e).  To  investigate  the  functional  role 
of  BRD4  hyperphosphorylation  in  BBI  resistance,  we  analysed  whether 
BRD4  phosphorylation  influences  MED1  binding.  Indeed,  SUM159R 
cells  treated  with  CK2  inhibitor  or  PTZ  both  lead  to  decreased  MED  1 
abundance  in  BRD4  immunoprecipitations,  suggesting  that  pBRD4 
binds  MED1  more  efficiently  than  BRD4  (Fig.  4f,  g). 

To  functionally  assess  the  role  for  BRD4  phosphorylation  in  BBI 
resistance  and  MED1  binding,  we  generated  BRD4  constructs  encoding 
mutants  that  cannot  be  phosphorylated  by  CK2  (7  serine  to  alanine 
substitutions;  £7A  mutant’)  or  mimic  constitutive  phosphorylation 
(7  serine  to  aspartate  substitutions;  £7D  mutant’).  We  first  assessed 
the  ability  of  these  constructs  to  rescue  effects  of  endogenous  BRD4 
knockdown  in  stable  cell  lines  (Extended  Data  Fig.  8b,  c).  We  observed 
expression  of  both  7D  and  7  A  mutants  supporting  the  growth  of  both 
parental  SUM  159  and  JQ1  -resistant  SUM159R  cells  (Fig.  3c).  Next, 
we  analysed  MED1  binding  and  subcellular  localization  of  7 A  and 
7D  mutants  in  the  presence  or  absence  of  JQ1.  We  found  that  the  7 A 
mutant  displays  weaker  MED1  binding  compared  to  wild-type  BRD4 
and  completely  dissociates  after  JQ1,  whereas  the  7D  mutant  seems 
to  have  higher  affinity  for  MED1  that  is  unaffected  by  JQ1  treatment 
(Fig.  4h  and  Extended  Data  Fig.  8i).  Lastly,  we  assessed  the  sensitivity 
of  cells  expressing  7  A  or  7D  mutant  BRD4  to  JQ1.  In  parental  SUM  159 
cells  exogenously  expressed  7D  mutant  BRD4  decreased  sensitivity  to 
JQ1,  whereas  the  7 A  mutant  slightly  increased  sensitivity  (Fig.  4i).  In 
contrast,  expression  of  7 A  mutant  BRD4  in  SUM159R  cells  restored 
JQ1  sensitivity  whereas  the  7D  mutant  showed  a  modest  decrease. 
These  results  strongly  support  the  hypothesis  that  hyperphosphoryl¬ 
ation  of  BRD4  arises  from  decreased  PP2A  activity  in  BBI  resistant 
cells,  leading  to  increased  binding  of  BRD4  to  MED1,  recruitment  to 
chromatin  and  decreased  responsiveness  to  bromodomain  inhibition. 

To  explore  the  clinical  relevance  of  phospho-BRD4  (pBRD4)  in  BET 
inhibitor-naive  TNBC,  we  performed  immunofluorescence  analysis  of 
a  tissue  microarray  (TMA)  featuring  of  89  patient-derived  TNBC  spec¬ 
imens.  First,  we  validated  the  pBRD4  immunofluorescence  assay  by 
comparing  xenografts  derived  from  SUM  159  and  SUM159R  cell  lines 
and  detected  significantly  higher  pBRD4  in  SUM159R  cells  (Extended 
Data  Fig.  9a).  We  detected  strong  pBRD4  staining  among  a  subset  of 
TNBCs  (Extended  Data  Fig.  9b,  c),  and  variable  staining  overall  that 
was  not  correlated  with  expression  of  the  androgen  receptor  and  basal 
cytokeratins  (bCK;  Extended  Data  Fig.  9d,  e)  and  it  was  not  signifi¬ 
cantly  associated  with  disease  outcomes  (Supplementary  Table  9  and 
Extended  Data  Fig.  9f). 

To  extend  the  translational  relevance  of  our  findings,  we  conducted 
synergy  studies  of  JQ1  with  molecules  targeting  BCL-xL  (ABT737),  a 
gained  super-enhancer  in  SUM159R  cells,  and  modulators  of  BRD4 
phosphorylation,  the  CK2  inhibitor  CX-4945  and  the  PP2A  activator 
perphenazine  (PPZ).  We  observed  significant  synergy  between  JQ1  and 
all  three  compounds  studied  (Extended  Data  Fig.  10),  establishing  a 
rationale  for  combination  studies  of  BBI  in  TNBC  to  improve  response 
and  to  anticipate  BBI  resistance. 

BRD4  inhibition  has  demonstrated  efficacy  in  disparate  models  of 
cancer  in  a  rapidly  expanding  literature.  Despite  apparent  resistance 
in  the  vast  majority  of  tumour  types,  as  we  observed  here  in  TNBC, 
mechanisms  of  BBI -resistance  have  not  been  mechanistically  explained. 
As  this  research  was  in  review,  two  studies  reported  moderate  emer¬ 
gent  resistance  to  BBI  in  murine  AML  associated  phenotypically  with 
a  stem-like  state  and  WNT  pathway  activation26,27.  Interestingly  in 
our  study  TNBCs  with  more  basal/ stem  cell-like  features  and  WNT 
pathway  activation  are  more  sensitive  to  BET  inhibition,  whereas 
resistant  disease  emerges  as  epigenomic  adaptation  to  a  more  differ¬ 
entiated  luminal  phenotype.  Our  findings  of  persistent  BET  bromodo¬ 
main  dependency  despite  BBI  resistance,  as  well  as  pBRD4  staining  in 
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resistant  disease,  should  be  studied  in  these  murine  AML  models  and 
further  in  human  leukaemia. 

Integrating  approaches  in  epigenomics,  proteomics  and  chemical 
biology,  we  provide  an  example  of  epigenomic  drug  resistance  by  an 
epigenetic  mechanism,  where  in  BBI -resistant  cells,  decreased  PP2A 
activity  leads  to  hyperphosphorylated  BRD4,  which  binds  more 
strongly  to  MED1,  facilitating  a  bromodomain-independent  chromatin 
recruitment  mechanism.  This  research  proposes  putative  combination 
strategies  to  anticipate  and  overcome  BBI  resistance,  including  pair¬ 
ing  with  BCL-xL  inhibitors  (for  example,  ABT-737)  or  CK2  inhibitors, 
and  guides  the  development  of  second-generation  BBIs  that  disrupt 
BET  function  via  orthogonal  biophysical  or  biochemical  actions.  More 
immediately,  the  robust  efficacy  observed  in  pre-clinical  models  sup¬ 
ports  the  development  of  BET  inhibition  in  TNBC  alone,  and  in  com¬ 
bination  with  mechanism-based  targeted  therapies. 

Online  Content  Methods,  along  with  any  additional  Extended  Data  display  items  and 
Source  Data,  are  available  in  the  online  version  of  the  paper;  references  unique  to 
these  sections  appear  only  in  the  online  paper. 
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No  statistical  methods  were  used  to  predetermine  sample  size.  The  experiments 
were  not  randomized,  and  the  investigators  were  not  blinded  to  allocation  during 
experiments  and  outcome  assessment. 

Cell  lines  and  breast  tumour  tissues.  Breast  cell  lines  were  obtained  from  the 
ATCC  and  S.  Ethier  (SUM  series).  Cells  were  cultured  in  media  recommended 
by  the  provider,  their  identity  was  confirmed  by  STR  analysis,  and  they  were  reg¬ 
ularly  tested  for  mycoplasma.  Breast  tumour  samples  were  collected  using  proto¬ 
cols  approved  by  the  DF/HCC  Institutional  Review  Board,  informed  consent  was 
obtained  from  all  patients.  Tumours  were  minced  with  razor  blades  and  digested 
with  stirring  for  3-4h  at  37°C  in  DMEM/F12  with  2  mg  ml-1  BSA,  2mgml_1  col- 
lagenase  type  IV,  and  2  mg  ml-1  hyaluronidase.  After  digestion,  cells  were  filtered 
through  500-pm  mesh,  washed  in  DMEM/F12  with  5%  FBS,  frozen  in  DMEM/ 
F12  with  5%  FBS  and  10%  DMSO,  and  stored  in  liquid  nitrogen  for  subsequent 
xenograft  studies.  PDX  IDC50  was  derived  from  a  primary  tumour  of  highly 
invasive  metaplastic  TNBC  resistant  to  chemo  and  radiation  therapy  leading  to 
the  rapid  death  of  the  patient.  Exome  sequencing  of  the  tumour  and  xenograft 
identified  numerous  mutations  including  heterozygous  frameshift  mutation  in 
PTEN  (chrl0_89701964-89701964_A)  and  CDH1  chrl6_67400242-67400242_C). 
PDX  EL- 12-58  was  derived  from  a  liver  metastasis  of  a  heavily  pretreated  basal-like 
TNBC,  Oncopanel  mutation  testing  identified  homozygous  mutations  in  BRCA2 
(p.S1970*),  TP53  (p.I232fs),  TSC2 ,  FLT3 ,  and  ROS1,  and  lower  frequency  muta¬ 
tions  in  RAD21,  JAK3,  ARID  IB,  ARID  1  A,  KDM6A. 

High-throughput  screening  of  BET  bromodomain  inhibitors  in  breast  cell  line 
panel.  We  tested  a  panel  of  compounds  (synthesized  in  the  Bradner  laboratory) 
in  40  human  breast  cell  lines  in  a  384-well  format  at  2,000  cells  per  well  using  a 
semi- automated  screen  essentially  as  described5.  Cell  viability  at  72  h  was  evaluated 
using  ATPlite  (Perkin  Elmer). 

Synergy  studies.  SUM149,  SUM149R,  SUM159,  and  SUM159R  cells  were  seeded 
in  sterile,  white,  opaque  384-well  microtitre  plates  (Thermo),  using  an  automated 
dispensing  system  (BioTek  EL406),  at  1,000  cells  per  well  in  50  pi  of  media.  Drugs 
were  delivered  in  DMSO  by  robotic  pin  transfer  with  a  JANUS  workstation  (100  nl) 
to  achieve  a  matrix  of  pairwise  dose-response  incubations  of  each  compound, 
each  pair  having  eight  replicates.  Following  72  h  of  incubation,  ATP  levels  were 
determined  for  treated  cells  and  vehicle  controls  (ATPlite,  PerkinElmer).  Data  were 
normalized  to  vehicle  controls.  Combination  indices  were  determined  using  the 
median-effect  principle  of  Chou  &  Talalay28  (CalcuSyn  Software).  Isobologram 
plots  were  generated  with  GraphPad  Prism  software.  Points  represent  paired  values 
of  drug  concentrations  assessed  for  synergism.  The  diagonal  line  signifies  drug 
additivity.  Points  above  the  line  represent  antagonistic  drug  combinations,  and 
those  below  the  line  represent  synergistic  drug  combinations.  Synergy  assays  were 
performed  in  triplicates  and  repeated  2-3  times. 

Xenograft  assays.  For  xenograft  assays  5-6-weeks  old  female  CrTacNCr -Foxnlnu 
and  NOD.Cg -Prkd(fcld  l72r/mJSwVJicTac  mice  were  purchased  from  Taconic. 
Tumours  were  induced  by  bilateral  orthotopic  mammary  fat  pad  injection  of 
1  x  106  cells  in  50%  Matrigel  (BD  Biosciences)  in  DMEM/F12  or  Medium  171 
(except  for  IDC50-X  cells,  which  were  injected  with  3%  FBS  and  4  mg  ml-1  col¬ 
lagen  gel  in  Medium  171).  Animal  experiments  were  conducted  following  proto¬ 
col  11-023  approved  by  the  Dana-Farber  Cancer  Institute  Animal  Care  and  Use 
Committee.  For  all  the  xenograft  studies,  the  sample  size  of  each  group  (5-10 
mice)  is  indicated  in  the  figures.  We  performed  pilot  experiments  using  a  few 
(5-10)  mice  per  group  followed  by  larger  studies  if  needed  to  reach  statistical 
significance  and  repeated  experiments  to  ensure  reproducibility.  Due  to  the  nature 
of  the  performed  experiments,  no  randomization  and  no  blinding  was  used  as  it 
was  deemed  unfeasible.  However,  the  resulting  tumours  were  analysed  in  a  blinded 
manner.  Mice  were  administered  JQ1  (50  mg  per  kg,  daily),  vehicle  only  (control) 
for  14  days  beginning  at  day  14  (SUM  159),  or  doxycycline  at  day  21  (SUM  159- 
shBRD4)  after  injection.  Mice  were  euthanized  and  tumours  evaluated  28  and 
60  days  after  injection  of  parental  and  TET-inducible  shBRD4-expressing  SUM  159 
cells  into  mammary  fat  pads. 

Cellular  viability,  senescence,  MDR  and  BH3  profiling  assays.  Cell  viability  and 
growth  assays  (Figs  la,  2a,  b,  3d,  e,  4c,  i,  Extended  Data  Figs  la,  b,  3i,  4d,  e,  g,  h,  10), 
cycle,  apoptosis,  and  MDR  assays  were  performed  in  triplicates  and  repeated  2-3 
times.  For  cell  proliferation  assays,  cells  were  plated  at  500  cells  per  well  in  96-well 
plates  and  treated  the  next  day  with  inhibitors,  DMSO  or  doxycycline  (500  ng). 
Cells  were  cultured  at  37  °C  with  5%  CO2  in  the  media  described  above,  and  cell 
viability  was  measured  using  CellTiter-Glo  three  days  after  treatments.  For  cell 
growth  assays,  cells  were  plated  at  5,000  (SUM159)  or  20,000  (SUM149)  cells  per 
well  in  6-well  plates  and  treated  the  next  day  with  inhibitors.  Cells  were  counted 
every  three  days  by  cell  counter.  Cellular  apoptosis  was  analysed  with  an  APC 
AnnexinV/7ADD  Apoptosis  Detection  kit  (BD  Pharmingen).  AnnexinV/7AAD 
assessments  and  cell  cycle  graphics  were  generated  using  Flowjo  software  V7.6.1 
for  Windows  (Tree  Star).  Senescence  Beta-gal  staining  was  performed  using 


Senescence  (3-Galactosidase  staining  kit  from  Cell  Signaling.  Briefly,  after  JQ1 
treatment  (500  nM)  for  72  h,  SUM159  and  MDA-MB-231  cells  were  fixed  by  fixa¬ 
tive  solution  for  15  min,  followed  by  p-galactosidase  solution  incubation  overnight 
at  37  °C.  The  staining  was  checked  under  microscope  for  the  development  of  blue 
colour.  Multi-Drug  Resistance  Assay  was  performed  with  MDR  assay  kit  from 
Cayman  Chemical  (600370).  Briefly,  SUM159  and  SUM159R  cells  were  treated 
with  JQ1  or  DMSO  for  30  min  in  SUM  medium.  Verapamil  was  used  as  a  positive 
control  at  1:1,000  dilution.  Calcein  AM/Hoechst  Dye  staining  solution  was  added 
after  that  and  cells  were  incubated  at  37  °C  for  15  min.  The  cells  were  analysed  by 
fluorescent  microscope  and  FACS. 

Cell  cycle  analysis  was  performed  72  h  after  JQ1  treatment  or  BRD4  down- 
regulation  with  doxycycline  using  propidium  iodide  (PI)  staining.  Cells  were 
resuspended  in  1  ml  of  growth  medium  supplemented  with  2pgml_1  PI  (Life 
Technologies)  as  final  concentration.  After  60  min  at  37  °C  in  the  dark,  analysis 
was  performed  on  a  FACS  Ariall  cytometer  (BD  Biosciences).  The  cell  cycle  was 
plotted  as  histogram  after  excluding  doublets. 

Cell  synchronization  procedure.  SUM  159  cells  were  treated  with  nocodazole 
(200  ng  ml-1)  for  12  h  and  then  cells  were  tapped  to  detach  from  the  plates.  After 
washing  twice  with  PBS,  cells  were  replated  with  or  without  JQ1  in  collagen-coated 
plates.  Cells  were  collected  at  0, 3, 6, 12  h  time  point  for  FACS  and  immunoblot  anal¬ 
ysis.  Dynamic  BH3  profiling  was  performed  using  the  JC-1  plate  method  as  previ¬ 
ously  described18,19.  Briefly,  2.5  x  105  cells  were  seeded  in  T25  flasks  in  the  presence 
of  500  nM  or  5  pM  JQ1  for  72  or  96  h.  Cells  were  trypsinized,  suspended  in  MEB 
(150mM  mannitol,  lOmM  HEPES,  50  mM  KC1,  5mM  succinate,  20  pM  EDTA, 
20  pM  EGTA,  0.1%  protease-free  BSA,  pH  7.5  ±0.1),  and  1-2  x  104  cells  were  added 
in  15  pi  of  MEB  to  each  well  of  a  384  well  Fluotrac  200  plate  containing  15  pi  per 
well  of  either  peptides  at  2  x  final  concentration,  buffer  only,  or  50  pM  alamethi- 
cin  in  MEB  supplemented  with  2  pM  JC-1,  lOmM  2-mercaptoethanol,  20pgml_1 
oligomycin,  and  50pgml_1  digitonin.  Fluorescence  at  emission  590  ±  10  nM  and 
excitation  545  ±  10  nM  was  recorded  at  5  min  intervals  at  32  °C.  The  area  under  each 
curve  was  normalized  to  the  alamethicin  and  buffer  controls  as: 

%  Depolarization  =1  [(AUC  sample  AUC  alamethicin)/ (AUC  buffer 
AUC  alamethicin)]  -  Delta  priming  was  calculated  per  peptide  treatment  as:  Delta 
Priming  =  (Depolarization  Treated)  —  (Depolarization  Untreated).  Positive  delta 
priming  indicates  an  increase  in  priming  due  to  treatment  and  an  increased  poten¬ 
tial  for  apoptosis  at  later  time  points. 

Immunofluorescence  staining  and  image  and  statistical  analysis  of 
tissue  microarrays.  Antibodies  used  for  immunofluorescence  were:  CK18 
(Dako,  M7010),  CK17  (Dako,  M7046),  HMW  CK  (Dako,  M0630),  LMW 
CK  (Dako,  M0631),  CD44  (NeoMarkers,  MS-668-P1),  CD24  (NeoMarkers, 
MS-1279-P1),  pSTAT3  (Cell  Signaling,  9145S),  VIM  (Dako,  M073501),  CDH1 
(BD  Biosciences,  610181),  Flag  (Sigma,  F1804),  BrdU  (Roche,  11170376001), 
pBRD4  (a  gift  from  C.  M.  Chiang),  and  androgen  receptor  (Cell  Signaling, 
5153S).  Immunofluorescence  experiments  were  performed  in  cultured  cells 
or  in  whole  sections  of  formalin-fixed  paraffin  embedded  (FFPE)  xenograft 
tumours.  The  staining  was  performed  as  described29.  Antibody  dilutions  were 
as  follows:  pSTAT3  (1:25),  CD44  (1:100),  CD24  (1:100),  CK18  (1:200),  CK17 
(1:200),  HMW  CK  (1:100),  LMW  CK  (1:100),  VIM  (1:100),  CDH1  (1:100), 
Flag  (1:50),  BrdU  (1:200),  pBRD4  (1:200),  and  androgen  receptor  (1:50).  The 
Dana-Farber  Breast  Cancer  Tissue  Microarray  (TMA)  consists  of  primary 
TNBC  samples  from  approximately  83  patients  who  underwent  definitive 
breast  surgery  at  Brigham  and  Womens  Hospital  between  1  January  1997  and 
31  December  2005.  Formalin-fixed,  paraffin-embedded  breast  cancers  were  col¬ 
lected  from  the  archives  of  the  Department  of  Pathology  at  Brigham  and  Womens 
Hospital.  Best  blocks  and  best  areas  for  coring  were  identified  and  selected  by  a 
breast  pathologist  (DD)  to  represent  different  area  of  the  tumour.  Results  of  immu- 
nohistochemical  studies  for  oestrogen  (ER)  and  progesterone  receptor  (PR)  and 
HER2  and  FISH  assay  results  for  HER2  were  extracted  from  pathology  reports. 
TMA  construction  was  carried  out  in  the  Dana-Farber/Harvard  Cancer  Center 
Tissue  Microarray  Core  Facility.  Three  0.6  mm  cores  were  taken  from  different 
marked  areas  in  each  case  and  placed  into  a  recipient  block  using  a  manual  arrayer 
(Beecher  Instruments).  Specimens  are  arrayed  in  triplicate.  Participants  signed 
consent  for  research  use  of  tissue  and  the  linking  of  tumour  specimens  to  clini¬ 
cal  follow-up.  Clinical  data  on  these  patients  was  collected  retrospectively  at  first 
presentation,  at  4,  9,  18,  30,  and  42  months,  and  annually  thereafter.  After 
9.3  years  median  follow  up,  24  recurrences  and  14  deaths  have  been  recorded. 
The  data  elements  are  the  following:  staging,  tumour  pathology,  diagnostic  and 
follow-up  tests  performed,  treatments  administered  (surgery,  radiation  and 
systemic  therapy),  and  recurrence.  Although  the  patients  in  this  cohort  were  not 
treated  as  part  of  a  clinical  trial  protocol,  they  were  treated  relatively  uniformly  as 
per  Dana-Farber  clinical  practice  guidelines.  This  serves  to  minimize  confound¬ 
ing  due  to  treatment  heterogeneity.  The  TMA  was  stained  with  pBRD4  (1:200) 
antibody  and  imaged  manually  on  Yokagawa  spinning  disc  confocal  microscope. 
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Three  images  were  taken  per  each  core  for  240  out  of 267  cores,  for  the  remaining 
27  one  or  two  images  were  taken  due  to  tissue  loss  or  low  tumour  content.  Image 
analysis  was  performed  with  Image J  software  macro  (code  available  upon  request). 
Phospho-BRD4  staining  mean  intensity  was  calculated  per  individual  nucleus 
within  an  image.  The  mean  intensity  per  image  was  normalized  to  nuclei  count. 
For  clinical  outcome  analysis  patients  were  dichotomized  as  ‘HighV’Low’  pBRD4 
by  median  intensity  (Supplementary  Table  9).  Disease  outcomes  were  evaluated  in 
83  of  89  TMA  samples  (3  were  not  TNBC  by  definitive  pathology,  2  did  not  have 
clinical  data  available,  one  was  a  repeat  biopsy  on  a  patient).  Recurrence-free  sur¬ 
vival  (RFS)  was  defined  as  the  interval  from  the  date  of  initial  surgical  resection  to 
the  date  of  recurrence  (local  or  distant),  or  date  of  last  known  contact  if  the  patient 
was  alive  and  has  not  recurred.  RFS  and  overall  survival  were  estimated  using 
the  Kaplan-Meier  product-limit  method,  with  hazard  ratios  and  95%  confidence 
intervals  from  a  univariate  Cox  proportional  hazard  model. 
siRNAs  and  lentiviral  shRNA  and  expression  constructs.  For  siRNA  transfec¬ 
tion  cells  were  plated  at  2,000  cells  per  well  in  96-well  plates  and  cultured  at  37  °C 
with  5%  C02  in  the  media.  The  next  day,  cells  were  transfected  in  triplicate  with 
siGENOME  SMARTpools  for  the  genes  of  interest  or  “Non-Targeting  siRNA”  con¬ 
trols  using  DharmaFECT  1  (Dharmacon).  The  sequences  of  the  siRNAs  in  the 
SMARTpools  are  listed  in  Supplementary  Table  10.  Cell  viability  was  measured 
using  CellTiter-Glo  (Promega)  three  days  after  transfections,  with  the  effects  of 
each  siRNAs  treatment  on  each  cell  line  compared  to  the  effects  of  no  siRNAs. 

TET-inducible  pLKO-TET-ON  lentiviral  constructs  were  packaged  by 
co-transfection  of  the  lentiviral  hairpin  containing  plasmid  PLKO.l  and  the 
helper  plasmids  pCMV-dR8.91  and  pMD2.G-VS.V-G  into  HEK293T  cells  using 
Lipofectamine  (Life  Technologies).  Following  transduction  via  spinoculation 
for  30 min  at  l,000g  and  selection  with  1  |igml-1  puromycin  for  72 h  (Sigma,  St. 
Louis,  MO),  knockdown  efficacy  was  determined  by  western  blotting  and  cells 
were  seeded  for  proliferation  assays  as  described  above.  Sequences  of  shRNAs  used 
are  listed  in  Supplementary  Table  10. 

Full  length  BRD4  in  pCDNA3  was  a  gift  from  Dr.  French  at  Brigham  and 
Womens  Hospital,  Harvard  Medical  School.  Mutations  of  BRD4  BD1  (N140A) 
and  BD2  (N433A)  bromodomains,  7A  and  7D  mutants  were  generated  using  a 
Quickchange  Multi  Site-Directed  Mutagenesis  Kit  (Agilent  Technologies)  using 
primers  listed  in  Supplementary  Table  10  and  subsequently  verified  by  sequencing. 
Immunoblotting  and  immunoprecipitation  experiments.  Cells  were  lysed 
five  days  after  transfection  with  siRNAs  in  RIPA  buffer.  Proteins  were  resolved 
in  SDS-polyacrylamide  gels  (4-12%)  and  transferred  to  PVDF  membranes  by 
using  a  Tris-glycine  buffer  system.  Membranes  were  blocked  with  5%  milk  pow¬ 
der  in  0.1%  Tween20  in  PBS  (PBS-T)  for  1  h  at  room  temperature  followed  by 
incubation  with  primary  antibodies  at  1:1,000  dilution  in  2.5%  milk  PBS-T.  For 
immunoprecipitation,  nuclear  extracts  were  prepared  as  follow:  1  x  107  cells  were 
resuspended  in  5  ml  buffer  A:  10  mM  Tris  pH  7.9, 1.5  mM  MgCl2,  10  mM  KC1, 
0.05%  NP-40,  1  mM  DTT,  and  protease  and  phosphatase  inhibitors.  Cells  were 
incubated  on  ice  for  15  min  and  gently  vortexed  every  5  min.  After  centrifugation 
at  2,000g  for  5  min,  pellets  were  suspended  in  0.3  ml  buffer  B  (20  mM  Tris  pH 
7.9,  25%  glycerol,  0.42  M  NaCl,  1.5  mM  MgCl2,  1  mM  KC1,  0.5%  NP40,  0.2  mM 
EDTA,  1  mM  DTT,  and  protease  and  phosphatase  inhibitors)  and  incubated  for 
5  min  on  ice.  After  centrifugation  of  the  lysates  at  14,000g  for  10  min  at  4  °C, 
supernatant  was  diluted  with  0.6  ml  buffer  A,  and  added  NP-40  to  final  0.5% 
and  treated  with  DNase  I.  The  samples  were  then  incubated  at  4°C  overnight 
with  BRD4  or  Flag  antibodies  at  1:100  dilution  and  immunoprecipitates  were 
collected  on  Dynabeads  Protein  G  for  2h.  Beads  were  washed  with  buffer  B 
containing  150mM  NaCl  and  0.5%  NP-40  three  times  and  then  resuspended 
in  gel  loading  buffer.  Immunoblotting  and  immunoprecipitation  experiments 
were  repeated  2-3  times. 

Antibodies  and  inhibitors.  Antibodies  used  for  immunoblotting,  immuno¬ 
precipitation  and  ChIP-seq  were  as  follows:  BRD4  (Bethyl,  A301-985A),  MED1 
(Bethyl,  A300-793a),  BRD3  (Bethyl,  A302-368A),  BRD2  (Bethyl,  A302-583A), 
MYC  (Santa  Cruz,  sc764),  pSTAT3  (Cell  Signaling,  9145S),  STAT3  (Cell  Signaling, 
4904),  pSTAT5  (Cell  Signaling,  9351),  pJAK2  (Cell  Signaling,  3771),  CYCLIN  D1 
(Cell  Signaling,  2922),  pH3  (Cell  Signaling,  12201),  CK2  substrate  (Cell  Signaling, 
8738),  PP2A-A  (Cell  Signaling,  2039),  PP2A-C  (Cell  Signaling,  2038)  and  pBRD4 
was  a  gift  from  C.  M.  Chiang.  Antibodies  used  for  ChIP-seq  were  BRD4  (Bethyl) 
histone  H3K27ac  (Abeam,  ab4729).  CXCR2  inhibitor  (239819)  and  CK2  inhibitor 
(218860)  were  from  CalBiochem,  JAK2  inhibitor  (Ruxolitinib,  INCB018424),  MEK 
inhibitor  (GSK1 120212,  S2673),  ABT-737  (sl002),  and  PI3K  inhibitor  (BKM120, 
S2247)  were  from  Selleckchem,  phenothiazine  (1525707)  and  perphenazine 
(1511000)  were  from  Sigma- Aldrich.  Inhibitor  treatment  for  immunoblot  analyses 
was  conducted  for  3  h. 

SILAC-RIME  experiments  and  data  analysis.  SUM159  and  SUM159R 
cells  were  grown  in  R/K-deficient  SILAC  DMEM  (paa;  E15-086),  10%  dialysed 
serum  (Sigma-Aldrich;  F0392),  and  supplemented  with  800  |iM  L-lysine  13C615N2 


hydrochloride  and  482  |iM  L-arginine  13C615N4  hydrochloride  (Cambridge  Isotope 
laboratory)  for  ‘heavy’ -labelled  media  or  800  |iM  L-lysine  12C614N2-hydrochloride 
and  482  pM  L-arginine  12C614N4  hydrochloride  for  ‘light’-labelled  media.  After 
SILAC  labelling,  RIME  was  performed  as  described20.  Word  clouds  for  Extended 
Data  Fig.  7  were  generated  using  R  version  3.1.0  and  the  R  package  ‘wordcloud’ 
version  2.5.  The  size  of  the  tag  reflects  the  square  root  of  the  MASCOT  score  of  the 
protein  (the  choice  of  square  root  is  arbitrary,  but  visually  appealing).  Experiments 
were  filtered  against  the  Contaminant  Repository  for  Affinity  Purification  Mass 
Spectrometry  Data30,  considering  any  protein  which  occurs  in  at  least  20  negative 
control  experiments  to  be  contamination,  hence  removed  from  the  data  set.  Refseq 
protein  IDs  provided  by  the  contaminant  repository  were  converted  to  Uniprot  IDs 
found  in  the  mass  spec  experiments  using  mappings  from  the  Bioconductor  pack¬ 
age  brg.Hs.eg.db’,  version  2.14  (Carlson  M.  org.Hs.eg.db:  Genome  wide  annotation 
for  Human.  R  package  version  3.0.0).  SILAC  RIME  experiments  were  performed 
in  duplicates  and  repeated  2-3  times. 

In  vitro  Chem-seq,  ChIP-seq,  and  RNA-seq.  Chem-seq  was  performed  essentially 
as  described11.  BRD4  ChIP-seq:  SUM159  and  SUM159R  cells  (4  x  107)  were  grown 
in  SUM  Medium.  The  media  were  then  removed  and  replaced  with  media  contain¬ 
ing  1%  formaldehyde  (EM  grade;  tebu-bio)  and  crosslinked  for  8  min.  Crosslinking 
was  quenched  by  adding  glycine  to  a  final  concentration  of  0.2  M.  The  cells  were 
washed  with  ice-cold  PBS,  harvested  in  PBS,  and  the  cell  pellet  was  washed  with 
PBS.  The  nuclear  fraction  was  extracted  by  first  resuspending  the  pellet  in  10  ml 
of  LB1  buffer  (50  mM  HEPES-KOH  (pH  7.5),  140  mM  NaCl,  ImM  EDTA,  10% 
glycerol,  0.5%  NP-40  or  Igepal  CA-630,  and  0.25%  Triton  X-100)  for  10  min  at 
4°C.  Cells  were  pelleted,  resuspended  in  10ml  of  LB2  buffer  (10 mM  Tris-HCL 
(pH  8.0),  200  mM  NaCl,  1  mM  EDTA,  and  0.5  mM  EGTA),  and  mixed  for  5  min. 
Cells  were  pelleted  and  resuspended  in  300  pil  of  LB 3  buffer  (lOmM  Tris-HCl 
(pH  8),  100  mM  NaCl,  1  mM  EDTA,  0.5  mM  EGTA,  0.1%  Na-deoxycholate,  and 
0.5%  N-lauroylsarcosine)  and  sonicated  in  a  Covaris  sonicator  for  10  min.  A  total 
of  30  pi  of  10%  Triton  X-100  was  added,  and  lysate  was  centrifuged  for  10  min  at 
20,000  ref  to  purify  the  debris.  The  supernatant  was  then  incubated  with  100  \A 
of  Dynabeads  Protein  G  (Life  Technologies,  10003D)  prebound  with  20  |ig  BRD4 
antibody  (Bethyl,  A301-985A),  and  immunoprecipitation  (IP)  was  conducted  over¬ 
night  in  the  cold  room.  The  beads  were  washed  ten  times  in  1  ml  of  RIPA  buffer 
and  twice  in  100  mM  ammonium  hydrogen  carbonate  (AMBIC)  solution.  DNA 
was  eluted  in  elution  buffer  (50 mM  Tris-HCl  pH  8, 10 mM  EDTA,  and  1%  SDS). 
Cross-links  were  reversed  overnight  at  65  °C.  RNA  and  protein  were  digested  with 
0.2  mg  ml-1  RNase  A  for  2  h  followed  by  0.2  mg  ml-1  Proteinase  K  for  1  h.  DNA  was 
purified  with  phenol- chloroform  extraction  and  ethanol  precipitation.  Libraries 
for  Illumina  sequencing  were  prepared  following  the  Rubicon  ThruPLEX-FD  kit 
for  10-12  cycles. 

RNA-seq:  SUM159  and  SUM159R  were  incubated  in  biological  duplicates  for 
3, 12  and  24  h  with  500  nM  of  JQ1  or  DMSO  treatment.  Total  RNA  was  extracted 
using  the  standard  Qiagen  RNeasy  kit  (74106).  RNA  concentrations  were  meas¬ 
ured  and  quality  controlled  on  a  Bioanalyzer,  RNA-Seq  libraries  were  made  using 
Illumina  True-Seq  RNA  kits  using  the  Sciclone  NGSx  workstation. 

All  RNA-seq  and  ChIP-seq  experiments  were  performed  in  duplicates. 
Genomic  data  analyses.  Accessing  data  generated  in  this  manuscript.  All  ChIP-seq, 
Chem-seq,  and  RNA-seq  data  generated  in  this  publication  can  be  found  online 
associated  with  GEO  Publication  Reference  ID  GSE63584  (www.ncbi.nlm.nih.gov/ 
geo/).  Supplementary  Table  2  lists  all  sequencing  data  sets  and  their  corresponding 
GEO  GSM  accession  IDs. 

Gene  sets  and  annotations.  All  analysis  was  performed  using  RefSeq  (NCBI37/ 
HG19)  human  gene  annotations. 

RNA-seq  data  processing  and  gene  expression  quantification.  All  RNA-Seq  data 
sets  were  aligned  to  the  transcriptome  using  Tophat231  (version  2.0.1 1)  using  the 
Illumina  igenomes  NCBI37/HG19  UCSC  transcriptome  build  retrieved  from 
http://ccb.jhu.edu/software/tophat/igenomes.shtml.  Alignments  were  performed 
using  default  parameters.  Transcript  expression  quantification  was  performed 
using  Cufflinks32  (version  2.2.0)  with  default  parameters  to  generate  gene  expres¬ 
sion  values  in  units  of  FPKM. 

ChIP-seq  and  Chem-seq  data  processing.  All  ChIP-seq  and  Chem-seq  data  sets 
were  aligned  using  Bowtie233  (version  2.2.1)  to  build  version  NCBI37/HG19  of  the 
human  genome  or  build  version  NCB37/MM9  of  the  mouse  genome.  Alignments 
were  performed  using  the  following  criteria:  -k  1,  with  all  other  parameters  set  to 
default.  These  criteria  preserved  only  reads  that  mapped  uniquely  to  the  genome 
without  any  mismatches. 

Calculating  read  density.  We  calculated  the  normalized  read  density  of  a  ChIP- 
seq  or  Chem-seq  data  set  in  any  region  using  the  Bamliquidator  (version 
0.9)  read  density  calculator  (https://github.com/BradnerLab/pipeline/wiki/ 
bamliquidator).  Briefly,  ChIP-Seq  reads  aligning  to  the  region  were  extended  by 
200  bp  and  the  density  of  reads  per  base  pair  (bp)  was  calculated.  The  density 
of  reads  in  each  region  was  normalized  to  the  total  number  of  million  mapped 
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reads  producing  read  density  in  units  of  reads  per  million  mapped  reads  per  bp 
(rpmper  bp). 

Identifying  ChIP-seq  and  Chem-seq  enriched  regions.  We  used  the  MACS  version 
1.4.2  (Model  based  analysis  of  ChIP-Seq)34  peak  finding  algorithm  to  identify 
regions  of  ChIP-Seq  enrichment  over  background.  A  P  value  threshold  of  enrich¬ 
ment  of  1  x  10-9  was  used  for  all  data  sets.  The  GEO  accession  number  and  back¬ 
ground  used  for  each  data  set  can  be  found  in  the  accompanying  Supplementary 
Table  2. 

Creating  heat  map  representations  of  ChIP-seq  occupancy.  Heat  maps  of  ChIP-seq 
occupancy  for  various  factors  were  created  as  described35.  Heat  maps  were  created 
for  the  ±10  kb  region  flanking  all  transcription  start  sites  (TSS)  or  for  the  ±10  kb 
region  flanking  all  TSS  distal  BET  bromodomain  bound  enhancers.  Each  row  plots 
a  specific  TSS  or  enhancer  region.  Rows  are  ranked  by  peak  occupancy  of  BET 
bromodomains  as  determined  by  Bio- JQ1  Chem-seq  signal  (Fig.  Id). 

Correlating  BRD4  and  H3K27ac  occupancy  to  Bio-JQl.  Occupancy  of  BRD4  and 
H3K27ac  was  correlated  to  Bio-JQl  occupancy  at  all  regions  of  Bio-JQl  enrich¬ 
ment  in  SUM  159  cells.  Pearson  correlation  statistics  are  shown  (Extended  Data 
Fig.  3a).  To  quantify  changes  in  BRD4  or  H3K27ac  occupancy  upon  JQ1  treat¬ 
ment,  all  Bio-JQl  enriched  regions  were  ranked  in  SUM159  cells  and  then  binned 
(n  =  10).  Corresponding  box  plots  of  BRD4  or  H3K27ac  log2  fold  change  with  or 
without  JQ1  are  shown  for  each  bin  (Extended  Data  Fig.  3b). 

Mapping  enhancers  and  super- enhancers  using  Bio-JQl  occupancy  or  BRD4. 
Enhancers  and  super  enhancers  (SEs)  were  mapped  using  the  ROSE  software  pack¬ 
age  described12,13  and  available  at  (http://younglab.wi.mit.edu/super_enhancer_ 
code.html).  In  SUM159  and  SUM159R  cells,  Bio-JQl  Chem-Seq  enriched  regions 
were  used  to  map  enhancers  and  super  enhancers  (Fig.  If).  In  SUM149  cells,  BRD4 
ChIP-seq  enriched  regions  were  used  to  map  enhancers  and  super  enhancers 
(Extended  Data  Fig.  3a).  Enhancers  are  defined  as  regions  of  Bio-JQl  binding  not 
contained  in  promoters. 

Quantifying  changes  in  gene  expression  of  super  enhancer  proximal  genes.  Genes 
within  50  kb  of  super  enhancer  in  SUM159  or  SUM149  were  identified  and  filtered 
for  expression  status  (>  1  FPKM  expression  in  any  sample),  and  filtered  to  remove 
non  poly-adenylated  transcripts  (for  example,  microRNAs).  For  SUM159,  log2  fold 
changes  in  gene  expression  at  super-enhancer-associated  genes  or  all  expressed 
genes  was  compared  at  3, 12,  and  24 h  post  JQ1  treatment  (Fig.  lg).  For  SUM149, 
comparisons  were  made  at  12  h  post  JQ1  treatment  (Extended  Data  Fig.  3f).  The 
statistical  significance  of  differences  between  distributions  of  changes  was  also 
assessed  using  a  Welch’s  two-tailed  t  test. 

Identifying  differentially  expressed  genes  upon  JQ1  treatment.  To  identify  genes 
differentially  regulated  by  JQ1  treatment  in  SUM  159  or  SUM  149  cells,  all  genes 
with  a  >  1  log2  fold  change  in  expression  were  ordered  by  fold  change  at  24  h  with 
or  without  JQ1  for  SUM159  or  at  12  h  with  or  without  JQ1  for  SUM149.  The  log2 
row  median  normalized  fold  change  for  each  gene  is  displayed  as  a  heat  map  in 
Extended  Data  Fig.  3g  for  SUM  159  and  in  Extended  Data  Fig.  3h  for  SUM  149. 
For  subsequent  gene  set  and  pathway  analysis,  SUM  159  genes  with  consistent  and 
statistically  significantly  altered  expression  were  selected  using  a  Welch’s  two-tailed 
t  test  between  DMSO  and  JQ1  treated  expression  values  at  12  and  24  h.  A  Rvalue 
cut-off  of  0.01  was  applied  (Extended  Data  Fig.  3k). 


Identifying  gained/lost  super  enhancers  between  SUM159  and  SUM159R.  Super 
enhancer  differential  regions  were  defined  as  in  Brown  et  al.  20 1436.  Briefly,  in 
order  to  quantify  changes  in  super-enhancers  between  two  conditions,  back¬ 
ground  subtracted  ChIP-Seq  signal  was  calculated  at  the  set  of  all  enhancer 
regions  considered  super  in  at  least  one  condition.  Gained/lost  super-enhancers 
were  determined  as  those  with  a  greater  than  log2  fold  change  signal  in  either 
direction.  The  log2  fold  change  in  Bio-JQl  occupancy  at  all  rank  ordered 
super-enhancer-containing  regions  is  shown  in  Fig.  2c.  Super  enhancer  regions 
were  classified  as  either  gained,  conserved,  or  lost.  Gained/lost  regions  were  clas¬ 
sified  as  those  with  >1  log2  fold  change  in  either  direction.  Conserved  regions 
were  classified  as  those  with  <0.25  log2  fold  change  in  either  direction.  The  log2 
fold  change  in  either  BRD4  or  proximal  (within  50  kb  of  region)  gene  expression 
is  shown  in  Extended  Data  Fig.  5d-f. 

Quantifying  changes  in  BRD4  and  H3K27ac  occupancy  upon  JQ1  treatment  in 
either  SUM  159  or  SUM159R  cells  at  Bio-JQl  regions.  Log2  fold  changes  in  BRD4 
or  H3K27ac  were  quantified  at  Bio-JQl  enriched  regions  in  their  respective  cell 
line  and  shown  in  Extended  Data  Fig.  5c. 

Quantifying  changes  in  BRD4  and  H3K27ac  as  a  function  of  Bio-JQl  or  BRD4 
occupancy.  Bio-JQl  enriched  regions  in  SUM159  or  BRD4  enriched  regions  in 
SUM  149  were  ranked  by  increasing  levels  and  then  distributed  into  10  bins.  Log2 
fold  changes  in  BRD4  or  H3K27ac  were  quantified  in  each  bin  of  regions  and 
displayed  as  a  box  plot  (Extended  Data  Fig.  3b,  c). 

Quantifying  changes  in  BRD4  occupancy  upon  JQ1  treatment  in  all  TNBC.  Log2  fold 
changes  in  BRD4  upon  JQ1  treatment  were  quantified  at  BRD4  enriched  regions 
in  each  respective  cell  line  Extended  Data  Fig.  6c. 

All  code  related  to  genomic  and  transcriptome  analysis  can  be  found  at  https:// 
github.com/BradnerLab/TNBC. 
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Extended  Data  Figure  1  |  See  next  page  for  figure  caption. 
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Extended  Data  Figure  1  |  BET  bromodomain  proteins  and  cell  growth 
in  TNBCs.  All  error  bars  represent  s.d.,  n  =  3.  a,  Cellular  viability  of 
SUM  159  and  MDA-MB-231  cells  expressing  TET-inducible  BRD4- 
targeting  or  lacZ  shRNAs.  P  values  indicate  statistical  significance  of  the 
observed  differences  (paired  t-test).  b,  Cellular  viability  four  days  after 
transfection  of  siRNAs  targeting  BET  bromodomain  proteins.  Asterisks 
indicate  statistical  significance  (paired  f-test)  of  the  marked  differences 
as  follows.  SUM159:  siBRD2  versus  siBRD3,  P  =  0.002;  siBRD3  versus 
siBRD4,  P  =  0.0006,  MDA-MB-231:  siBRD2  versus  siBRD3,  P  =  0.006; 
siBRD2  versus  siBRD4,  P  =  0.002;  siBRD3  versus  siBRD4,  P  =  0.016, 
MDA-MB-468:  siBRD2  versus  siBRD3,  P=  0.0009;  siBRD3  versus 
siBRD4,  P  =  0.0055,  MDA-MB-436:  siBRD2  versus  siBRD4,  P=  0.002; 
siBRD3  versus  siBRD4,  P=  0.015,  ZR-75-1:  siBRD2  versus  siBRD3, 

P=  0.0169;  siBRD3  versus  siBRD4,  P  =  0.007.  c,  Immunoblot  analysis  of 
BET  bromodomain  proteins  four  days  after  siRNA  transfection. 


d,  Cell  cycle  profile  of  SUM  159  cells  synchronized  in  G2/M  with 
100  ng  ml-1  nocodazole  followed  by  replating  to  fresh  medium  with 
DMSO  or  JQ1  (500  nM)  added  at  —  1  h  or  at  3  h  after  release.  Cells  were 
collected  at  different  time  points  (0,  6,  12  h)  after  release,  e,  Immunoblot 
analysis  of  the  indicated  proteins  at  different  time  points  (0,  3,  6,  12  h) 
after  release  of  SUM159  cells  synchronized  with  lOOngml-1  nocodazole 
followed  by  replating  to  fresh  medium  with  DMSO  or  JQ1  (500  nM)  added 
at  1  h  before  or  3  h  after  release,  f,  Cell  cycle  analysis  of  SUM  159  cells 
following  72  h  treatment  with  JQ1  (500  nM)  or  downregulation  of  BRD4 
using  TET-inducible  shRNAs.  g,  Annexin  V  staining  of  SUM  159  cells 
following  72  h  treatment  with  JQ1  (500  nM)  downregulation  of 
BRD4  using  TET-inducible  shRNAs.  All  error  bars  represent  s.e.m. 
h,  Immunoblot  analysis  of  the  indicated  proteins  in  a  panel  of  breast  cell 
lines;  colour  scheme  as  in  a.  For  gel  source  data,  see  Supplementary  Fig.  1. 
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Extended  Data  Figure  2  |  See  next  page  for  figure  caption. 
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Extended  Data  Figure  2  |  Response  to  BBIs  in  TNBCs.  a,  Immunoblot 
analysis  of  the  indicated  proteins  at  different  time  points  following 
JQ1  treatment  (500  nM)  in  SUM159  cells  (top)  and  at  different  JQ1 
doses  for  24 h  treatment  in  SUM159  and  MDA-MB-436  cells  (bottom), 
b,  Immunoblot  analysis  of  the  indicated  proteins  at  different  time 
points  following  JQ1  treatment  (500  nM)  in  SUM149,  SUM159  and 
MDA-MB-231  cells,  c,  Haematoxylin  and  eosin  staining  of  SUM159  cells 
after  3  days  of  JQ1  treatment,  d,  Senescence  (3-galactosidase  staining  of 
SUM159  and  MDA-MB-231  cells  after  3  days  of  JQ1  treatment.  Scale 
bars,  100  [im.  e,  Box  plots  depict  the  weights  of  xenografts  30  days  after 
injection  of  MDA-MB-231  (2  x  106)  and  IDC50X  (2  x  105)  cells  into 
inguinal  mammary  fat  pads  of  NOG  mice;  n  indicates  the  number  of  mice 
per  experiment.  P  values  indicate  statistical  significance  of  the  observed 
differences  (unpaired  £-test).  Error  bars  represent  s.e.m.  Mice  were 
administered  JQ1  (50  mg  per  kg,  daily)  or  vehicle  only  (control)  for  14 
days  beginning  at  day  16  (MDA-MB-231)  or  10  (IDC50X)  after  injection 
(after  tumours  reached  palpable  size).  For  EL12-58X  PDX,  mice  were 


implanted  with  pieces  of  tissue  measuring  1x3x3  mm  into  the  inguinal 
mammary  fat  pads  and  were  administered  daily  JQ1  (50  mg  per  kg)  for 
14  days  beginning  at  day  21  after  injection  (after  tumours  reached  palpable 
size),  f,  Bromodeoxyuridine  (BrdU)  and  luminal  (Low  MW  CK)  and  basal 
high  (High  MW  CK)  molecular  weight  cytokeratin  staining  of  EL12-58 
xenograft  with  or  without  JQ1  treatment.  Scale  bars,  50  [im.  g,  Tumour 
volume  of  SUM  159  cells  expressing  TET-inducible  BRD4-targeting 
shRNAs.  Mice  were  administered  doxycycline  or  vehicle  only  (control) 
for  39  days  beginning  at  day  21  after  injection  (after  tumours  reached 
palpable  size).  Error  bars  represent  s.d.,  n  =  4  (shBRD4-l  experiment) 
and  n  =  5  (shBRD4-2  experiment),  h,  Haematoxylin  and  eosin  staining 
and  immunofluorescence  analysis  of  basal  (basal  cytokeratin,  cytokeratin 

17,  pSTAT3,  and  CD44)  and  luminal  (luminal  cytokeratin,  cytokeratin 

18,  and  CD24)  markers  in  SUM159  xenografts  with  or  without  JQ1 
treatment.  Scale  bars,  100  \im  for  haematoxylin  and  eosin  and  50  \im  for 
immunofluorescence,  respectively.  For  gel  source  data,  see  Supplementary 
Fig.  1. 
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Extended  Data  Figure  3  |  SUM149  JQ1  response,  a,  Scatter  plots  showing 
the  relationship  between  the  genomic  binding  of  BRD4  and  Bio-JQl 
(left)  or  H3K27ac  and  Bio-JQl  (right)  at  all  Bio-JQl  enriched  bound 
regions.  Units  of  genomic  occupancy  are  in  rpm  per  bp.  A  simple  linear 
regression  is  drawn  in  black.  Pearson  correlation  statistics  are  also  shown, 
b,  Box  plots  showing  the  log2  fold  change  in  BRD4  with  or  without  JQ1 
(left)  or  H3K27ac  with  or  without  JQ1  (right)  at  Bio-JQl  bound  regions 
in  SUM159  cells.  The  12,999  Bio-JQl  regions  are  ranked  by  increasing 
Bio-JQl  binding  and  divided  into  10  separate  bins  (displayed  from  left 
to  right).  The  statistical  significance  of  the  difference  in  the  mean  BRD4 
log2  fold  change  between  the  weakest  and  strongest  Bio-JQl  bound  region 
bins  is  shown  (Welch’s  t-test;  ***P  value  <1  x  10-10).  c,  Box  plots  showing 
the  log2  fold  change  in  BRD4  with  or  without  JQ1  (left)  or  H3K27ac 
with  or  without  JQ1  (right)  at  BRD4  bound  regions  in  SUM149  cells. 

The  5,696  BRD4  bound  regions  are  ranked  by  increasing  background 
subtracted  BRD4  binding  and  divided  into  10  separate  bins  (displayed 
from  left  to  right).  The  statistical  significance  of  the  difference  in  the 
mean  BRD4  log2  fold  change  between  the  weakest  and  strongest  BRD4 
bound  region  bins  is  shown  (Welch’s  f-test;  ***P  value  <1  x  10-10).  d, 
Ranked  plots  of  enhancers  defined  in  untreated  SUM  149  cells  ranked  by 
increasing  BRD4  signal  (units  rpm).  Enhancers  are  defined  as  regions  of 
BRD4  binding  not  contained  in  promoters.  The  cut-off  discriminating 
typical  from  super-enhancers  is  shown  as  a  dashed  grey  line.  Enhancers 
associated  with  TNBC  characteristic  genes  are  highlighted,  e,  Scatter 
plots  showing  the  relationship  between  the  log2  fold  change  in  gene 
expression  upon  12  h  JQ1  treatment  in  SUM149  (y-axis)  and  SUM159 
(x-axis)  cells.  A  simple  linear  regression  is  drawn  in  red.  The  Pearson 


correlation  statistic  is  also  shown,  f,  Box  plots  showing  the  log2  fold  change 
in  expression  relative  to  DMSO  control  of  either  all  active  genes  or  super¬ 
enhancer  (SE)-associated  genes  upon  12  h  JQ1  treatment.  The  statistical 
significance  of  the  difference  in  expression  change  between  all  active 
genes  and  super-enhancer-associated  genes  is  shown  by  a  Welch’s  t-test  *P 
value  <  1  x  10-3).  g,  Heat  map  showing  the  expression  of  genes  that  are 
up  or  down  regulated  by  JQ1  versus  DMSO  after  24  h  treatment.  Each  row 
shows  the  expression  of  a  single  gene  in  either  DMSO  or  JQ1  treated  cells 
at  3,  12,  and  24  h  after  treatment.  Expression  values  are  coloured  according 
to  fold  change  relative  to  the  median  for  each  row.  Genes  are  ordered  by 
fold  change  with  or  without  JQ1  24  h  after  treatment,  h,  Heat  map  showing 
the  expression  of  genes  that  are  up  or  down  regulated  by  JQ1  versus  DMSO 
after  12  h  treatment  in  SUM149  and  SUM149R  cells.  Each  row  shows  the 
expression  of  a  single  gene  in  either  DMSO  or  JQ1  treated  cells  at  12  h  after 
treatment.  Expression  values  are  coloured  according  to  fold  change  relative 
to  the  median  for  each  row.  Genes  are  ordered  by  fold  change  with  or 
without  JQ1  12  h  after  treatment  in  SUM149  cells,  i,  j,  Box  plots  showing 
the  log2  fold  change  in  expression  at  genes  that  are  up-  (i)  or  down- 
(j)  regulated  by  JQ1  versus  DMSO  after  12  h  of  treatment  in  parental 
SUM149  cells.  Log2  fold  change  in  expression  is  shown  for  either  parental 
SUM149  (left)  or  resistant  SUM149R  (right)  cells,  k,  Top  signalling 
pathways  affected  by  JQl-induced  gene  expression  changes  in  SUM159 
cells.  1,  Viable  cell  numbers  of  SUM  149  (left)  and  SUM149R  (right) 
cells  treated  with  different  doses  of  JQ1  (2  [xM,  10  fxM).  Error  bars 
represent  s.d.,  n  =  3.P  values  indicate  statistical  significance  of  the 
observed  differences  (two-way  ANOVA  with  Bonferroni  multiple 
comparison  correction). 
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Extended  Data  Figure  4  |  See  next  page  for  figure  caption. 
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Extended  Data  Figure  4  |  Characterization  of  SUM159R  cells. 

a,  Expression  of  ABC  transporters  in  SUM159  and  SUM159R  cells.  The 
expression  of  29  ABC  transporters  was  analysed  based  on  RNA-seq  data 
on  the  two  cell  lines,  b,  Assay  for  MDR  (multi  drug  resistance)  pumps 
in  SUM159  and  SUM159R  cells  treated  with  JQ1  alone  or  together  with 
verapamil  based  on  microscopic  examination  (left)  and  FACS  (right)  of 
cells  labelled  with  fluorescent  MDR  substrate,  c,  Immunoprecipitation 
analysis  of  Biotinylated  JQ1  (Bio-JQl)  in  SUM159  and  SUM159R  cells 
with  JQ1  treatment  at  different  time  points  following  immunoblot  for 
the  indicated  proteins,  d,  Cellular  viability  of  SUM  159  and  SUM159R 
cells  treated  with  CXCR2  and  JAK2  inhibitors.  Error  bars  represent  s.d., 
n  —  3.  e,  Cellular  viability  of  SUM  159,  and  pool  and  single  cell  clones  of 
SUM159R  cells  treated  with  different  doses  of  JQ1.  Error  bars  represent 
s.d.,  n  =  3.  f,  Tumour  weight  of  xenografts  derived  from  SUM  159  and 


SUM159R  cells.  Mice  were  administered  JQ1  for  14  (SUM159)  and  30 
(SUM159R)  days  beginning  at  day  14  and  26,  respectively,  after  injection. 
P  values  indicate  statistical  significance  of  the  observed  differences 
(unpaired  t- test).  Error  bars  represent  s.e.m.  g,  Immunoblot  analysis 
of  BCL-XL  expression  in  SUM159  and  SUM159R  cells  before  and  after 
JQ1  3  h  treatment  (500  nM).  h,  Dynamic  BH3  profiling  reveals  inverse 
apoptotic  response  to  JQ1  in  SUM149R  and  SUM159R  cells.  In  parental 
lines  JQ1  increases  priming  relative  to  untreated  cells  indicating  an 
increase  in  apoptotic  propensity.  In  resistant  lines  JQ1  reduces  priming 
indicating  greater  resistance  to  apoptosis  relative  to  untreated  cells. 

P  values  indicate  statistical  significance  of  the  observed  differences 
(two-way  ANOVA).  Error  bars  represent  s.e.m.,  n  —  5.  For  gel  source  data, 
see  Supplementary  Fig.  1. 
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Extended  Data  Figure  5  |  See  next  page  for  figure  caption. 
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Extended  Data  Figure  5  |  BRD4  binding  in  SUM159R  cells,  a,  Cellular 
viability  of  SUM  159  and  SUM159R  cells  transfected  with  siRNAs  targeting 
bromodomain  proteins.  Asterisks  indicate  statistical  significance  (paired 
t-test)  of  the  marked  differences  as  follows:  SUM  159:  siBRD2  versus 
siBRD3,  P  =  0.013,  siBRD3  versus  siBRD4,  P  =  0.0154  and  SUM159R: 
siBRD2  versus  siBRD3,  P  =  0.0159,  siBRD2  versus  siBRD4,  P  —  0.0048; 
siBRD3  versus  siBRD4,  P=  0.0068.  b,  Cellular  viability  of  SUM159R  cells 
expressing  TET-inducible  BRD4-targeting  or  lacZ  shRNAs.  All  error  bars 
represent  s.e.m.  Rvalues  indicate  statistical  significance  of  the  observed 
differences  (unpaired  f-test).  c,  Box  plot  showing  the  log2  fold  change  in 
H3K27ac  genomic  occupancy  at  regions  bound  by  Bio-JQl  in  parental 
SUM  159  or  resistant  SUM159R  cells,  d,  Heat  map  showing  the  expression 
of  genes  that  are  up  or  down  regulated  by  JQ1  versus  DMSO  after  24  h 
treatment  in  parental  SUM  159  cells.  Each  row  shows  the  expression  of  a 
single  gene  in  either  DMSO  or  JQ1  treated  cells  at  24  h  after  treatment  in 
SUM  159  cells  (left  four  columns)  or  SUM159R  cells  (right  four  columns). 
Expression  values  are  coloured  according  to  fold  change  relative  to  the 
median  for  each  row.  Genes  are  ordered  by  fold  change  with  or  without 
JQ1  24  h  after  treatment,  e,  f,  Box  plots  showing  the  log2  fold  change 
in  expression  at  genes  that  are  up-  (e)  or  down-  (f)  regulated  by  JQ1 


versus  DMSO  after  24  h  of  treatment  in  parental  SUM  159  cells.  Log2  fold 
change  in  expression  is  shown  for  either  parental  SUM  159  or  resistant 
SUM159R  cells,  g,  Box  plots  showing  expression  of  genes  that  are  up-  or 
downregulated  by  JQ1  versus  DMSO  after  24  h  of  treatment  in  parental 
SUM  159  cells.  Expression  is  shown  in  DMSO  and  JQ1 -treated  conditions 
in  units  of  FPKM  for  either  parental  SUM  159  (left)  or  resistant  SUM159R 
(right)  cells.  The  statistical  significance  of  the  difference  between  gene 
expression  distributions  for  SUM159  DMSO-  and  JQl-treated  cells  is 
shown  (P  <  0.01).  The  difference  between  all  other  distributions  are 
considered  non  significant  (NS).  The  statistical  significance  of  the 
difference  between  SUM  159  DMSO  gene  expression  distribution  and 
all  other  distributions  is  shown  (*P  <  1  x  10-3).  The  difference  between 
all  other  distributions  are  considered  non-significant,  h,  Examples  of 
luminal  and  basal  cell- specific  genes,  and  MYC  in  SUM  159  and  SUM159R 
cells.  RNA-seq  tracks  are  shown,  i,  Haematoxylin  and  eosin  staining  and 
immunofluorescence  analysis  of  luminal  (CK18  and  LMW)  and  basal 
(CK17  and  HMW)  cytokeratins  and  luminal  (VIM  and  CD24)  and  basal 
(CDH1,  CD44,  and  pSTAT3)  cell  markers  in  SUM159R  xenografts.  All 
error  bars  represent  s.e.m.  Scale  bars,  100  \im  for  haematoxylin  and  eosin 
and  50  \xm  for  immunofluorescence,  respectively. 
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Extended  Data  Figure  6  |  JQ1  response  in  other  breast  cancer  cell  lines. 

a,  b,  Gene  tracks  depicting  BRD4  +  DMSO  and  BRD4  +  JQ1  in  multiple 
breast  cancer  cells  at  the  BCL-xL  (a)  or  SOD2  (b)  gene  loci.  The  x-axis 
shows  position  along  the  chromosome  with  gene  structures  drawn  below. 
They-axis  shows  genomic  occupancy  in  units  of  rpm  per  bp.  The  BCL-xL 
and  SOD2  super- enhancers  are  shown  as  a  red  bar  at  the  top.  c,  Box  plots 
showing  the  log2  fold  change  in  BRD4  occupancy  with  or  without  JQ1  for 


all  BRD4  bound  regions  in  each  cell  line  for  multiple  TNBC.  Cell  lines 
are  ordered  by  their  JQ1  (IC5o)  and  coloured  by  their  sensitivity,  d,  Gene 
tracks  depicting  H3K27AC  occupancy  at  the  BCL-xL  locus  in  SUM149 
parental  (top,  light  blue)  or  SUM149R  resistant  (bottom,  dark  blue)  cells. 
The  x-axis  shows  position  along  the  chromosome  with  gene  structures 
drawn  below.  They-axis  shows  genomic  occupancy  in  units  of  rpm  per  bp. 
All  error  bars  represent  s.e.m. 
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Extended  Data  Figure  7  |  Word  clouds  depicting  BRD4-associated  proteins  identified  in  RIME  analysis. 
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Extended  Data  Figure  8  |  Mechanism  of  BBI-resistance.  a,  Immunoblot 
analysis  of  BRD4  immunoprecipitates  for  MED1  in  the  indicated  cell  lines 
with  or  without  JQ1  treatment  (5  pM,  3  h).  b,  Immunoblot  analysis  of  long 
(BRD4L)  and  short  (BRD4S)  forms  of  BRD4  after  transfection  of  siRNAs. 
c,  Immunoblot  analysis  of  the  indicated  exogenously  expressed  Flag- 
tagged  BRD4  proteins  in  SUM  159  and  SUM159R  cells,  d,  Immunoblot 
analysis  of  pBRD4,  BRD4,  MED1  and  ACTB  in  the  indicated  cell  lines 
with  or  without  JQ1  treatment,  e,  Immunoblot  analysis  of  phospho-BRD4 
(pBRD4)  and  BRD4  in  SUM  159  and  SUM159R  cells  treated  with  the 
indicated  doses  of  CK2,  PI3K,  and  MEK  inhibitors  for  2  h. 


f,  Immunoblot  analysis  of  CK2  substrates  in  SUM  159  and  SUM159R  cells 
following  CK2  inhibitor  (CX-4945,  10  pM)  3  h  treatment,  g,  Immunoblot 
analysis  of  pBRD4  and  BRD4  in  SUM  149  cell  line  treated  with  different 
doses  of  the  indicated  PP2A  inhibitors  for  3  h.  ACTB  was  used  as 
loading  control,  h,  Immunoblot  analysis  of  pBRD4  and  BRD4  in  the 
indicated  cell  lines  treated  with  different  doses  of  phenothiazine  for  6h. 
i,  Immunofluorescence  analysis  of  exogenous  Flag-tagged  BRD4  proteins 
(WT,  BD,  7D  and  7A)  in  SUM159  cells  with  or  without  JQ1  treatment 
(5  pM,  3  h).  Scale  bars,  20  pm.  For  gel  source  data,  see  Supplementary 
Fig.  1. 
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Extended  Data  Figure  9  |  See  next  page  for  figure  caption. 
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Extended  Data  Figure  9  |  Phospho-BRD4  levels  in  xenografts  and 
primary  TNBC  samples,  a,  Immunofluorescence  analysis  of  phospho- 
BRD4  (pBRD4)  in  SUM159  parental  and  SUM  159R  xenografts  showing 
that  resistance  is  associated  with  higher  pBRD4  levels,  b,  Examples  of 
pBRD4  immunofluorescence  in  patient  tumours  depicting  variability 
among  different  TNBC  samples.  Scale  bars,  50  [im.  c.  Mean  intensity  of 
phospho-BRD4  (pBRD4)  in  tissue  samples  from  83  patients  with 
early- stage  triple  negative  breast  cancer  (TNBC).  d,  Examples  of  androgen 
receptor  and  basal  cytokeratin  (bCK,  HMW  CK)  immunofluorescence 
in  TNBC  samples.  Scale  bars,  50  \im.  e,  Box  plot  depicting  pBRD4  signal 
intensity  in  TNBCs  tumours  with  the  indicated  androgen  receptor  and 
bCK  expression  patterns.  None  of  the  differences  among  groups  were 


significant  (ANOVA  test  P  =  0.5413  and  Dunnetf  s  multiple  comparisons 
test  was  not  significant),  f,  Kaplan-Meier  estimates  of  disease-free  survival 
(DFS)  and  overall  survival  (OS)  in  TNBC  subgroups  using  a  median-split 
of  pBRD4  intensity.  Disease  outcomes  were  evaluated  in  83  of  89  TMA 
samples.  Patients  with  low  pBRD4  had  a  worse  overall  prognosis  with 
a  five-year  RFS  of  66.2%  (95%  confidence  interval  (Cl)  52.7-83.1%), 
compared  to  an  RFS  of  86.4%  (95%  Cl  76.0-98.3%)  among  patients  with 
high  pBRD4  (HR  =  2.3,  95%  Cl  0.98-5.4,  P  =  0.06).  However,  with  this 
small  sample  size  this  difference  did  not  reach  statistical  significance,  nor 
did  the  association  of  pBRD4  status  and  overall  survival  (HR  =  2.0,  95%  Cl 
0.67-5.9,  P  =  0.22). 
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Extended  Data  Figure  10  |  Overcoming  BBI-resistance.  a-c,  Synergy  The  diagonal  line  signifies  drug  additivity.  Points  above  the  line  represent 

studies  of  JQ1  with  ABT737  (BCL-xl  and  BCL-2  inhibitor)  (a),  CX-  antagonistic  drug  combinations,  and  those  below  the  line  represent 

4945  (CK2  inhibitor)  (b)  and  perphenazine  (PP2A  activator)  (c).  Points  synergistic  drug  combinations, 
represent  paired  values  of  drug  concentrations  assessed  for  synergism. 
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